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1.0  INTRODUCTION 


Over  the  last  17  years  many  uses  of  TMA  (trimethyl 
aluminum)  trails  in  the  upper  atmosphere  have  been  made  by 
various  experimenters  (Blamont  et  al,  1961;  Armstrong,  1963; 
Rosenberg  et  al,  1964).  In  addition  to  measuring  winds  and 
diffusion  rates,  (Authier  et  al,  19G2,  1963,  1964;  Authier, 
1964;  Harang,  1964;  Golomb  et  al,  1967,  1972;  Fees  et  al, 

1974)  have  also  measured  temperature  by  studying  the  tempera¬ 
ture  dependent  vibrational-rotational  intensity  distribution 
of  the  blue-green  system  of  A10  (B2Z+  -  X2Z+)  generated  when 
TMA  is  exposed  to  sunlight  in  the  presence  of  atomic  oxygen 
(see  Figure  1). 

Recent  experiments  under  the  direction  of  AFGL 
(Aladdin,  Aeolus,  etc.)  have  centered  around  rocket  released 
puffs  of  TMA  on  the  downleg  of  a  rocket  trajectory  from  189  km 
down  to  120  km  with  approximately  10  km  spacing  (see  Figure  2) 

A  program  involving  one  rocket  payload  for  winds  and 
temperature  studies  was  conducted  at  White  Sands  Missile  Range 
on  24  September,  1977.  In  this  case  the  rocket  made  only  one 
release  at  approximately  160  km  altitude  at  5:00  A.M.  local 
time  during  morning  twilight  (8°  solar  depression).  The  ex¬ 
perimental  procedure  consisted  of  observing  the  top  and  bottom 
of  the  cloud  for  30  to  40  second  periods.  In  general,  ade¬ 
quate  signal  to  noise  ratio  is  maintained  until  1  or  2°  solar 
depression  although,  after  this  length  of  time  (approximately 
20  minutes)  from  an  initial  release  at  say  8°  solar  depression 
only  at  the  lowest  altitudes  120-130  km  are  there  "sharp" 
enough  features  to  permit  altitude  discrimination  in  the  spec¬ 
trometer's  line  of  sight  and  accurate  photographic  triangula¬ 
tion  of  altitude.  In  this  program,  only  12  minutes  of  data 
was  taken  with  the  above  limitation  being  reached  more  quickly 
at  the  higher  altitudes  available  in  this  particular 
experiment . 
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2.0  INSTRUMENTATION  AND  PROCEDURE 


The  heart  of  field  instrumentation  is  a  1/2  meter 
Fastie-Ebert  scanning  spectrometer.  The  input  optics  is  a 
50  cm  focal  length  F/5.0  lens  which  "images”  the  input  slit 
upon  the  sky  as  a  horizontal  slit  7°  wide  by  .06°  high.  The 
spectrometer  is  co-mounted  on  a  tripod  adjustable  in  eleva¬ 
tion  and  azimuth  along  with  a  secondary  electron  conduction 
(SEC)  vidicon  low  light  level  television  camera  (and  a  35  mm 
photographic  camera  as  backup)  for  both  boresight  recording 
and  real-time  aiming  of  the  spectrometer.  The  grating  of  the 
spectrometer  is  tilted  to  scan  in  wavelength  by  a  synchronous 
motor  through  an  adjustable  cam  mechanism  at  a  rate  of  one 
scan  every  five  seconds.  Photoelectron  pulses  from  the  photo¬ 
multiplier  are  processed  with  standard  "nuclear"  counting 
electronics  and  recorded  on  a  high  speed  (60  ips)  magnetic 
tape  recorder.  Also  recorded  are  a  sample  of  the  60  Hz  wave¬ 
form  driving  the  synchronous  motor  scanning  the  spectrometer 
and  audio  comments  of  the  person  aiming  the  spectrometer.  A 
block  diagram  of  the  Field  Instrumentation  is  shown  in  Figure  3. 

The  boresight  TV  camera  is  recorded  on  a  video  tape 
recorder  along  with  time  from  a  digital  clock  superimposed  on 
one  corner  of  the  video  frame.  Also  recorded  on  the  video 
tape  recorder's  audio  channel  is  the  commentary  of  the  "aimer" 
describing  what  he  is  pointing  the  instrument  at  as  well  as 
when  he  is  moving  it  or  holding  steady  on  some  feature  (see 
Figure  4). 

The  aimer  makes  use  of  a  TV  monitor  for  aiming  in 
addition  to  various  direct  sights  mounted  on  the  tripod.  He 
also  has  view  of  a  chart  recorder  displaying  the  spectrum  be¬ 
ing  recorded  in  order  to  ascertain  the  quality  (amplitude) 
of  the  spectra  from  a  particular  feature. 
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Bores ight  Cameras  I Video 


?  Operator 
Lc rc phone 


figure  M  . 


Bores ight  TV  Monitor  (Typical) 


After  the  experiment  certain  calibrations  are  re¬ 
corded.  An  iron-neon  hollow  cathode  lamp  is  scanned  to  es¬ 
tablish  the  wavelength  calibration  of  the  spectrometer  and  the 
shape  of  the  spectrometer's  slit  function.  A  light  box  (LDF 
#6)  with  a  radiance  calibration  traceable  to  an  NBS  quartz- 
iodine  irradiance  standard  is  scanned  to  establish  the  varia¬ 
tion  of  the  spectrophotometer's  sensitivity  with  wavelength. 
Records  of  the  photomultiplier  output,  with  the  spectrometer 
scan  stopped,  versus  the  step  number  of  the  calibrated  aper¬ 
tures  in  the  lightbox  are  used  to  establish  the  amplitude 
linearity  of  the  counting-recording  system.  In  a  photon  count¬ 
ing  scheme  such  as  this,  counts  are  lost  at  higher  counting 
rates  due  to  "dead-time".  Prior  to  the  experiment  boresight 
calibrations  using  a  several  kilometer  distant  light  source 
are  made  to  align  all  the  sights  as  well  as  markers  on  the 
boresight  TV  monitor  and  fiduciaries  on  the  photographic 
camera  with  the  very  narrow  (vertically)  field  of  view  of  the 
spectrometer.  Later,  during  data  reduction,  playback  of  the 
TV  recording  of  this  boresight  calibration  will  again  allow 
us  to  place  markers  on  the  TV  monitor  used  for  playback. 


3 . 0  DATA  REDUCTION 


Data  reduction  was  carried  out  in  the  following 
steps.  A  block  diagram  of  data  reduction  is  shown  in  Figure  5. 

1.  Digitize  Data  -  The  tapes  brought  back  from  the  field 
are  "replayed”  and  the  photon  counts  are  totaled  every 
l/60th  of  a  second  using  as  a  clock  the  recorded  60  Hz 
waveform  originally  driving  the  synchronous  motor  on  the 
grating  drive.  The  totalized  counts  are  recorded  on  a 
digital  tape  recorder  whose  operation  as  well  as  the 
counting  and  interface  to  the  analog  "playback"  tape 
recorder  is  handled  by  a  microprocessor  based  system. 

The  following  programs  are  used  on  the  AFGL  CDC  6600  computer: 

2.  Program  CONVERS  -  The  "stranger"  type  digital  tape 
created  above  is  converted  to  a  SCOPE  format  file  by 
CONVERS. 

3.  Program  ALO  TAPE  DUMP  -  Prints  all  of  the  data  on  the 
high  speed  printer  in  an  easy  to  examine  format  together 
with  a  running  index  to  keep  track  of  position  along  the 
tape  (see  Table  1). 

4.  Program  LB  AVG  takes  sections  of  data,  that  have  been 
identified  by  examining  the  dump  above  as  observations 
of  the  lightbox  at  various  aperture  steps,  and  computes 
the  average  count  corresponding  to  the  lightbox  steps  0, 

1,  2,  4,  8,  16,  .  .  .  .,  1024  (see  Table  2). 

5.  Program  LBTAU  takes  the  averages  computed  above  (k1) 
plus  the  measured  "hole  factors"  for  each  attenautor 
step  on  the  light  box  (k)  and  adjusts  the  parameters 

a,  0,  t  in  the  expression  below  for  best  fit  to  properly 
weighted  values  of  k  and  k' . 
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where  k'  =  measured  counts 

k  =  output  of  lightbox 
x  =  dead  time 

a  =  constant  of  proportionality  since  the  hole 
factors  for  the  aperture  steps  of  the  light¬ 
box  are  in  arbitrary  units,  not  counts 

0  -  an  additive  constant  representing  the  dark 
current  of  the  photomultiplier. 

From  this  point  forward  in  the  data  reduction  the 
data  can  be  corrected  to  remove  the  non-linearity  introduced 
by  the  dead  time  inherent  in  the  counting  by  the  expression 


where  k"  is  the  corrected  counts  (see  Figure  6). 

6.  Program  KITQPL  makes  plots  on  the  CRT  plotter  of  the  raw 
data  once  the  starting  points  of  the  scans  are  established 
from  examining  the  dump  of  data  in  ALO  TAPE  DUMP.  These 
"quick  plots"  can  be  visually  correlated  with  the  video 
tape  boresight  record  by  comparing  the  relative  amplitude 
of  the  quick  plot  with  the  appai’ent  brightness  of  the 

A10  image  in  the  field  of  view.  Once  this  timing  has 
been  established,  a  table  is  created  of  sequential  scar, 
numbers  on  the  data  tape  versus  the  identification  num¬ 
bers  on  the  video  boresight  record  (see  Figure  7). 

7.  Program  KITSTPL  makes  strip  plots  on  computer  printout 

of  the  delta  tapes.  Groups  of  sequential  scans  of  the  same 
feature  in  the  TMA  release  are  summed.  Summed  scans  of 
the  lightbox  are  plotted  this  way  to  establish  the  varia¬ 
tion  of  sensitivity  wilh  wavelength  of  tie  instrument 
using  the  known  calibration  curves  of  the  lightbex 
(traceable  to  KBS).  This  slight  variation  ol  sensitivity 


M  ( COUNTS/SEC  ) 


M  COUNTS/SEC ) 


Figure  6.  Dead  Time  Curves 
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(~5%  over  100  A)  is  used  later  to  correct  the  data  in 
MATCH  (see  Figure  8). 

Scans  of  the  iron  hollow  cathode  lamp  lines  establish 
the  wavelength  calibration  of  the  spectrometer  system 
as  well  as  the  actual  shape  of  the  slit  function. 

Program  KITSLIT  takes  the  summed  scans  across  selected 
lines  and  fits  a  triangular  slit  function  to  them,  adjust¬ 
ing  an  arbitrary  multiplicative  constant,  a  background 
pedestal,  the  horizontal  position  and  the  width.  Then 
scans  of  different  lines  can  all  be  normalized  and  plotted 
on  computer  printout  as  shown  in  Figure  9 .  An  accurate 
slit  function  can  be  deduced  from  this  composite  by  hand 
drawing  smoothed  sections  at  the  top  and  bottom  of  the 
straight  line  of  the  triangular  function  (physically,  the 
sharp  corners  of  the  geometric  optically  correct  triangu¬ 
lar  slit  function  are  rounded  off  by  diffraction). 

Program  SYN  calculates  a  synthetic  A10  spectra.  The 
technique  is  similar  to  Harang  (1966).  The  molecular 
parameters  of  Tyte  and  Nicholls  (1964)  are  used  to  cal¬ 
culate  the  position  of  the  rotational-vibrational  lines. 
Populations  of  various  rotational  and  vibrational  levels 
are  derived  from  the  Boltzmann  distribution  law  assuming 
thermal  equilibrium  with  the  ambient  air  which  would  ap¬ 
pear  to  be  a  valid  assumption  up  to  approximately  170  km 
where  the  collision  rate  is  10/sec;  the  f-value  of  the 
A10  electronic  transition  is  generally  assumed  to  be  0.1 
to  1.0  (Harang,  1966)  implying  that  each  molecule  has  1  to 
10  seconds  or  10  to  100  collisions  to  attain  thermal 
equilibrium  between  excitations.  The  solar  excitation 
used  is  taken  from  the  Utrecht  Atlas  (Minnaert  et  al, 
1940).  The  hand  strength  with  Franck  Condon  factors  is 
used  to  calculate  absorption  for  various  vibrational 
levels . 
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Assuming  dynamic  equilibrium  between  vibrational  levels 
and  between  rotational  levels  relative  intensities  can 
be  assigned  to  each  rotational  line  in  the  bands  of  a 
sequence . 

The  intensities  of  the  first  75  rotational  lines  of  each 
branch  (P  and  R  branches)  in  the  first  three  bands  of  the 
AV  =  0  sequence  are  then  convolved  with  the  experimentally 
measured  slit  function  derived  in  KTTSLIT.  This  convolu- 

O 

tion  is  then  tabulated  every  1/4  A  and  every  50  degrees 
from  150°K  to  1050°K  and  stored  as  a  permanent  file. 

Plots  of  synthetic  spectra  are  shown  in  Appendix  B. 

.  Program  MATCH  is  the  heart  of  the  data  reduction  scheme. 
Inputs  to  MATCH  are  synthetic  spectra  from  SYN  already 
convolved  with  the  slit  function  at  19  temperatures  and 
the  digitized  spectrometer  output.  Groups  of  scans  are 
summed  or  lumped  as  determined  by  the  examination  of 
quick  plots  from  KITQPL  and  the  study  of  the  boresight 
records.  The  spectrometer  output  is  then  corrected  for 
amplitude  non-linearity  or  dead-time  with  subroutine 
FKAY  using  constants  supplied  by  program  KITAU.  Correction 
is  then  made  for  variation  of  sensitivity  with  respect  to 
wavelength  using  a  linear  approximation  derived  from  the 
calibration  scans  of  the'  lightbox. 

The  data  points  in  the  lumped  measured  scans  are  then 
compared  to  each  of  the  19  synthetic  spectra.  Iterative 
adjustments  are  made  to  the  exact  horizontal  (wavelength) 
position  of  the  measured  scan,  a  multiplicative  factor 
which  scales  the  synthetic  spectra  to  the  measured  scan, 
and  an  additive  constant  which  is  added  under  the  syn¬ 
thetic  spectra  as  a  pedestal.  The  iterative  adjustments 
of  the  three  variables  described  above  and  shown  graphi¬ 
cally  in  Figure  10  are  made  using  as  a  criterion  for  best 
fit  of  the  measured  points  to  the  synthetic  curve  the 
likelihood  that  the  measured  points  result  from  a  curve 
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MEASURED  SPECTRUM 


ADJUSTMENTS  TO  OBTAIN  OPTIMUM  FIT 
«  ADDITIVE  CONSTANT- BACKGROUND  PEDESTAL 

-  MULTIPLICATIVE  CONSTANT -SCALING  FACTOR 

-  HORIZONTAL  SHIFT- WAVELENGTH  ADJUSTMENT 


Figure  10.  Adjustments  Used  in  MATCH 
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same  as  the  synthetic  but  made  noisy  according  to  Poisson 
statistics.  The  likelihood  function  is  the  product  of 
the  probabilities  at  each  point  in  the  measured  scan 
that  the  point  results  from  Poisson  distributed  errors 
about  a  point  at  the  same  wavelength  in  the  synthetic 
spectrum.  Incidentally,  this  type  of  maximum  likelihood 
function  procedure  simplifies  to  become  a  least  square 
fit  when  the  errors  in  the  measurement  are  Gaussian.  It 
has  been  determined  that  in  this  application,  involving 
counting  randomly  arriving  photons  which  is  properly 
described  by  Poisson  statistics,  with  a  fairly  low  number 
of  counts  at  each  point  in  the  measured  spectrum,  con¬ 
sidering  that  Gaussian  statistics  become  a  poorer  approxi¬ 
mation  to  Poisson  statistics  at  smaller  counts  with  a 
consistent  bias  from  the  skewed  nature  of  the  Poisson  dis 
tribution  and  that  the  temperature  deduced  by  this  tech¬ 
nique  varies  monotonically  with  the  slope  of  the  data 
points,  use  of  least  squares  fitting  rather  than  the  maxi 
mum  likelihood  technique  based  on  Poisson  statistics  re¬ 
sults  in  a  systematic  error  of  as  much  as  15  degrees. 

The  temperature  of  the  synthetic  spectrum  producing  the 
maximum  likelihood  function  as  well  as  the  next  higher 
and  next  lower  temperature  and  the  three  corresponding 
values  of  the  likelihood  function  are  used  to  generate  a 
parabolic  curve  for  a  parabolic  interpolation  (see  Figure 
11).  The  peak  of  this  parabola  is  used  as  the  interpo¬ 
lated  temperature  of  best  fit.  More  details  are  supplied 
in  Appendix  A. 

In  Appendix  B  plots  of  the  fourteen  "lumped"  or  summed 
scans  of  measured  spectra  are  shown  along  with  five  syn¬ 
thetic  spectra  closest  in  temperature  to  each  of  the 
measured  spectra.  The  measured  spectra  are  corrected 
for  dead-time  amplitude  non-linearity  and  the  variation 
of  instrumental  sensitivity  with  respect  to  wavelength. 
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where: 

N-NO.  OF  INDEPENDENT  MEASUREMENTS 

STD.  DEVIATION  OF  TEMPERATURE  MEASUREMENT 
Ls  SUM  OF  SQUARES  (OR  EQUIVALENT  FOR  POISSON 
STATISTICS) 

T* TEMPERATURE 


Figure  11.  Interpolation  of  MATCH  Results  (Example) 


The  measured  spectra  as  well  as  the  five  synthetic 
spectra  are  normalized  to  peak  at  100  on  the  two  decade 
1  to  100  scale.  The  synthetic  spectra  are  shown  "sit¬ 
ting  on  top  of"  a  pedestal  C.j/k^  where  and  are  the 
constants  derived  by  the  fitting  procedure  for  best 
match  of  the  particular  synthetic  spectrum  to  the  mea¬ 
sured  spectrum.  Even  though  only  the  region  from  4824 
o 

to  4886  A  is  used  in  the  fitting  procedure,  more  of  the 
measured  spectrum  including  part  of  the  spectrometer  re¬ 
turn  sweep  is  shown  on  the  plot. 
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4.0  RESULTS 


The  results  of  the  September  1977  White  Sands  temper¬ 
ature  measurements  are  summarized  in  Table  3.  Since  there  was 
only  one  puff  or  cloud,  the  data  was  taken  from  points  near 
the  top  and  near  the  bottom  of  the  cloud  rather  than  from  the 
center  of  each  in  a  series  of  small  puffs  as  would  usually  be 
done  in  this  sort  of  experiment.  An  average  temperature  for 
the  observations  of  the  top  of  the  cloud  weighted  by  l/o2  is 
622°K  with  an  overall  o  of  22°K.  Disregarding  the  three  weak 
measurements  where  there  may  have  been  instrumental  problems, 
the  average  is  628°K  with  a  o  of  16°K.  An  average  temperature 
for  the  observations  of  the  bottom  of  the  cloud  with  or  with¬ 
out  the  two  weak  measurements  was  578°K  with  an  overall  o  of 
14°K. 

Some  of  the  scatter  of  individual  measurements, 
particularly  at  the  bottom  of  the  cloud,  is  a  result  of  having 
in  the  line  of  sight  to  the  bottom,  the  more  rapidly  diffusing 
upper  part  of  the  cloud.  In  general,  only  the  earliest  mea¬ 
surements  of  A10  at  altitudes  as  high  as  these  are  accurate 
for  reasons  of  altitude  spread  along  the  line  of  sight  through 
the  cloud. 

Better  altitude  resolution  could  have  been  had  if 
the  elevation  angle  of  the  line  of  sight  were  to  have  been 
lower;  however,  in  the  case  of  the  White  Sands  site  the  logistic 
requirements  precluded  moving  the  observation  site  back  100- 
200  km  as  would  have  been  preferred.  Also  the  accident  of 
having  a  very  large  release  over  perhaps  7-10  km  rather  than 
the  smaller  "point"  releases  did  not  aid  the  altitude  preci¬ 
sion  either. 
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WSMR  A10  TEMPERATURE 
24  September  1977 


File  No.  2  AMUL  =  0.3010 


Calculations  19  June  78 
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APPENDIX  A 


STATISTICAL  DEVELOPMENT  OF  FITTING  PROCEDURE 


APPENDIX  A 


The  general  likelihood  function  where  the  error 
distributions  at  each  point  are  independent  is 


Li 


n 

n 

i  =  l 


Aij 


i'Vi> 


(1) 


where  =  table  of  values  representing  measurements  at 

successive  points  in  the  measured  spectrum. 


P^ijCyi>  =  the  probability  that  a  point  in  the  spectrum 
having  a  long  term  average  A^  would  be  mea¬ 
sured  as  y^  in  a  single  measurement. 


A 


ij 


k.  z. .  +  C. 
1  11  3 


(2) 


where  z ^  =  table  of  values  of  synthetic  spectra  at  wavelengths 
corresponding  to  y^  at  the  same  i  index  and  tem¬ 
perature  j . 

kj  =  an  arbitrary  multiplicative  factor  used  to  scale 
the  synthetic  spectrum  to  the  measured  spectrum 
and  will  be  optimized  to  obtain  the  maximum  like¬ 
lihood  at  a  given  j . 

Cj  =  an  arbitrary  additive  constant  added  to  the  syn¬ 
thetic  spectrum  to  represent  a  constant  background 
or  pedestal  under  the  measured  spectrum.  It  will 
be  optimized  same  as  k^ . 

In  this  case  where  the  distribution  of  individual 
measurements  about  a  long  term  average  is  a  Poisson  distribution 


Aij 
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Combining  (1),  (2)  and  (3)  and  taking  the  logarithm  of  both 


n 

tn  Lj  =  .J  [y.  4n(kj  z +  C  j  )  -  (kj  z^  *  C . )  -  in(y!)]  (4 ) 

It  is  necessary  now  to  find  the  maximum  likelihood  of  in  Lj 
as  a  function  of  k  j ,  Cj  and  j  independently.  The  third  term 
in  (4)  is  not  a  function  of  the  above  so  it  can  be  ignored 
for  the  time  being.  At  any  given  j  then  the  optimum  kj  and 
Cj  are  those  which  allow  the  partial  derivatives  below  (5), 
(6)  to  be  zero 
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Unfortunately  there  is  no  analytic  solution  for  kj  and  Cj  as 
in  the  analogous  procedure  for  a  Gaussian  distribution  (least 
squares).  The  Newton  Rapsohn  method  will  be  used  to  effect  a 
solution. 
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in  general 
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solving  the  right-hand  side  of  both  equations. 
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new  values  for  k 

and  C  are 
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(15) 
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+  6C 

(16) 

which  are  iterated  until  6ic/k  and  6C/C  are  both  less  than 
10"  a  or  so.  Convergence  to  the  10"#  criteria  takes  five  or 
six  iterations. 


Once  optimum  kj  and  values  are  found  correspond¬ 
ing  to  a  particular  the  actual  likelihood  function  at  each 
j  can  be  computed.  The  at  which  the  likelihood  function 
is  maximum  is  the  temperature  of  best  fit.  The  preceding  T ^ 
along  with  T.+^  and  ^j  +  2  are  use<^  to  generate  a  parabola 


L  =  AT2  +  BT  +  C  (17) 

The  peak  of  the  parabola,  where  3L/3T  =  0  is  used  to  indicate 
an  interpolated  value  of  T.  The  shape  of  the  parabola  near 
the  peak  is  used  tc  deduce  a  in  a  manner  similar  to  Dyer 
(1970,  unpublished)  for  least  squares  fitting  (see  Figure  11). 
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APPENDIX  B 

VISUAL  VERIFICATION  OF  MATCH  RESULTS 
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. . . .  i  nnwrffr 


rn 

M 


xiiSN'nNi  '■jMiHVjy 


.Q  70  62  32  4  900 

F.NC'tH  (  AnOSt<:'OmO  i 

SEQUENCE  <hQ 


0  y  i  s  o  n  $ ;  h  .i'jci  j  1 3  a  y  m 


SCAN  NOS. 


SCAN  NOS  •  2 


UlSN'-UNl  IAiih'3M 

*2 


SCON  NOS. 


UISN'UNI  3AIiH13y 


PROGRAM  ASP < OU T  PUT , TAPES , TAPCA , TAPE3, TAPE 7 r TAPE 39 ) 

CC  PROGRAM  TMEO 

DIMENSION  TE<  2) »W(2) *UX(2) » WY( 2> *iE(2>  * ALFA( 2) . I Z ( 20  > , I IZ ( 20  > 
DIMENSION  FC<8*0> .R(0*8)*8NV(0) * SNJ< 0*80) »EMI (0*00) *APS<0*0O> 
DIMENSION  FR ( 3*3*00 >,FP< 3* 3, 00), JEX(0> 

DIMENSION  MJR <  00 > .HJ0(80> *HJP(0O> *0LAP<0O) , BIP<80) *BLAR(BO) * 
IDIR(BO) 

COMMON  TE  ,  W,  WX  »WY  *  BE  *  ALFA 
COMMON/ SPEC TP/ SLAP • PL AR »BIP*BIR 
WRITE  <  6*  29/ ) 

297  FORMAT  (38H1AL0  PROGRAM  i960  FOLLOWING  DATA  USED  > 

REWIND  7 

READ  <  3 . 293 ) 1  DATE » MOG 
UR  I TE ( A »  293 1 ) I DATE* MOG 

293  FORMAT  (16,12) 

2951  FORMAT  <1X*I6»X2) 

READ  (5,296)<JEX(IJEX),I JEX- 1*3) 

WRITE<6,296>( JEX(UEX)* IJ£X-1*5> 

294  FORMAT < 1H  »SA10> 

READ ( 5*109) TSOL  * YSOL  * SOLO 
WRITE(6»1091 ) TSOL  , TSOL » SOLO 
109  FORMAT (2F10.0*E12.4) 

1091  FGRMAT(1X*2F10.0*E12.4) 

R£AD<5»111 )(TE(I) , U< I ) ,UX< I > » UY i I >*BE(I) *  ALFA< I > » 1-1 ,2 > 

UR I TE(6,1111)(TE(I)»U(I>  »UX< I),WY(I),BE(I),  ALFA( I ) * 1*1 * 2) 

111  FORMAT! F12.4,F11 . 5 *F1 2 . 7  * F9 . 6 , F9 . 7 , FI  0 . 7 > 

1 11 1  FORMAT (lX»F12*4fF11.3»F12»7rF9.4tF9.7fF10»7) 

RE  AD ( 5 » 1 1 2  >  0  *  I MAX » JM AX 
WR!TE<6* 1 121 )G» IMAX* JMAX 

112  FORMAT (F20. 3*213) 

1121  F0RMAT(1X*F20. 3*213) 

DO  113  !•! t IMAX 

DO  113  1 1  *  1  *  I MAX 

REAP(5*114)L*LL»FC(I*II),R(I*II> 

WRITE! 6, 1141 )L»LL*FC(I,II>*R( J.II) 

114  FORMAT  (2I2»E1 1 • 4*F7 » 3 ) 

1141  F0RMAT(1X*2I2»E11.4,F7.3> 

113  CONTINUE 

READ<5* 122>H* 1AX , XX  AX*  JNAX 
UR  I rE( 4* 1221 )M# IAX » I IAX r JNAX 

122  FORMAT (413) 

1221  FORMAT  ( 1 X  *  413 ) 

IF(M>127, 127*1233 

127  DO  128  I-lrlAX 
DO  128  II-1,IIAX 
DO  128  J-l  ► JNAX 

REAP! 5* 1 80) IA » IB*LQ*  JP  *  TP# FP( I , 1 1  *  J) ,  JO, JR  *  TR » FR( I » 1 1 »  J ) 

WRITE  <6* 1001 > IA» 1 B  »LG*  JP, TP»FP<I*IIrJ) *  JO, JR » TR*FR< I  *  X  X  •  J ) 

100  FORMAT (2I3*I4*I3*F10t2*F10.3*2l3,F10.2*F10.3> 

1801  FORMAT<lX»2I3»I4tI3.F10.2rF10*3r2I3»F10.2fFl0.3) 

128  CONTINUE 

1233  READ ( 5» 193>LL  * JKAX 

WRITE(6*1?31 )LL*  JKAX 

READ  ( 5 » 193  )(IZ!L>»IIZ!L)*L-1»LL) 

WRI TE(4»193l)<IZ<L)» 1 1 7  <  L ) »L-1 *LL  ) 

193  FORMAT  (1414) 

1931  FORMAT (1X*14I4) 

READ! 5*1 15 >T. CELT. TMAX 
WRITE (4*1151 ) T »  DEL  T »  TMAX 
TBEGIN=T 

115  FORMAT  (3F10.1) 

1151  FORMAT (IX* 3F 10*1) 

292  FORMAT  < 24H  THESE  DATA  ON  TAPE  IF  1  -,I2> 

WRITE  <  4*  292) MOG 

READ! 5 *  1 93  ) ISPECP » IL I NE *  ICON, l INV 
WRITE (4* 1931 ) ISPECP* IL 1NE  *  ICON* 1 1 NO 

123  CONTINUE 

GT= 1 . 0/<0.6952f T ) 

JAX=*  JMAX  +  1 
V-0.0 
0-1.0 

114  V-VM.O 

00- ( W( 2 ) -WX <  2 > )*U-WX!2)*U*V 
EX*EXF  (  GOUOT) 

0-01EX 

IF (EX* 200  <0-1. 0)117*117.11* 

117  OV-O 

DO  120  L • t  *  I MAX 
VL-L  -  1 

GOO* ( W( 2  > -WX ( 2 ) >  *VL -UX  <  2  >  *Ut  *VL 
UN-( 1 ,0/0V)*rxr  ( -GOOtfiT ) 

SNVU  )*VN 
120  CONTINUE 

R-DFC*>  -Al  rAC’jtl  .3 


000720 

000730 

000740 

000750 

000760 

000770 

000780 

000790 

000800 

000010 

000820 

000830 

000840 

000850 

000860 

000870 

000880 

000090 

000900 

000910 

000920 

000930 

000940 

000950 

000960 

000970 

000900 

000990 

001000 

001010 

001020 

001030 

001040 

001050 

001060 

001070 

001080 

001090 

001100 

001110 

001120 

001130 

001140 

001150 

001160 

001170 

001180 

001190 

001200 

001210 

001220 

001230 

001240 

001250 

001260 

001270 

001280 

001290 

001300 

001310 

001320 

001330 

001340 

001350 

001360 

001370 

001380 

001390 

001400 

001410 

001420 

001430 

001440 

001450 

001460 

001470 

001400 

001490 

001500 

001310 

001520 

001530 

001540 

001550 


D*4.Ot0F(2>t»C(2>lK<2>/<IJ(2>tW<Z>  > 

AJ-0.0 

Q**  1  .0 

146  AJ-AJ41 . 0 
A-AJt<AJ+1.0> 

FJ-BM-D*A*A 
F-EXP  ( -F J*GT ) 

FA-<AJ4AJ41 .0>*F 
0-Q4F  A 

JF<  < FA* 100.0/0) -1 .0)147*147.146 

147  OJ-Q 

DO  148  I  *  1  *  I MAX 
V-I-l 

BB-BE  <2>-ALF A<2>*<04  0.5> 

DO  140  J-l.JAX 
SJ-J'l 

S- ( S  JF 1 . 0  >  4SJ 
FS-BP*S-P*S*S 
F-EXP  ( -FStOT ) 

SF-(SJ+SJ41 .OJIF/DJ 

148  SN  J  < I  *  J ) -SF 
B0-0.69524TS0L 
SOL -SLOF  <  YSOL  » SO ) 

C-SOLO/SOL 

DO  353  1= 1 » I MAX 
DO  150  J-l » JHAX 
SUH=O,0 

DO  152  IJ-l#IhAX 
UR-1 • O-GAR  < I # 1 1 ) 

FK-UR*UR*FC< I » 1 1  > 

SU«SNV<IJ> 

SIY-SU4FK 

TYR-SNY (I»II»J»J-1 > 

TYP«SNY< I » 1 1 » J» J+&  > 

ZF(M)179f 179» 154 
179  IF< I-IAX)174*l74tl54 
174  IF < II-IIAX) 170* 178» 134 
170  IF<  J- JHAX) 176* 176 » 154 
176  FRR-FR<t*IIr J) 

FRP»FP< I » 1 1  *  J ) 

GO  TO  153 
154  FRR= 1 « 0 
FRP=1 .0 

133  IF< J-l > 155t 153* 156 
153  SIR-0.0 
00  TO  157 

156  SJ«SNJ<II#J-l> 

CJ-J-t 

CCJ-2*< J-l )+l 
SR-CJ/CCJ 
5L-C6SL0F ( TYR.SO) 
SIR«SJ*SR*TYR*SL*FRR 

157  SSJ-SNJUIrJft) 

DJ-J 

DDJ-24J41 

SP-BJ/DDJ 

SSL=C4SL0F (TYP» SO) 

SIP-S5J*SP*SSL*TYP*FRP 

SI-SIU*(SIP+SXR> 

SUN-SUH+SI 
152  CONTINUE 

APS(If J)-SUH 
150  CONTINUE 
355  CONTINUE 

00  160  1-1 » 7NAX 
PO  161  J-l t JHAX 
SUM-0.0 

DO  162  1 1-l » INAX 
UR-1  .0*GIR(MI) 

FK*URfURtFC( I » II ) 

IF ( J-l )  163. 163*164 

163  SIR-0.0 
00  TO  165 

164  SR- J-l 

SYR-8NY  < I » 1 1  *  J  J* 1 ) 

S3R*  SYRtSYRAfcYP 
SIR-SR6S3R 

165  SP-J 

SYP-SNY  < I » 1 1 f  J»  J41 ) 

S3P-SYP4SYP4SYP 
SIF-SPBS3P 
SI*<Sir'4SIR»4FK 
162  SUH-SUM4SI 
161  ENT ( I  * J)«SUM 


001560 
001570 
001380 
001590 
001600 
001610 
001620 
001630 
001640 
001650 
001660 
001670 
001600 
001690 
0017 00 
001710 
001720 
001730 
001740 
001750 
001760 
001770 
001780 
001790 
00 1000 
001B10 
001820 
001830 
001840 
001850 
001860 
001870 
001880 
001090 
001900 
001910 
001920 
001930 
001940 
001950 
001960 
001970 
001980 
001990 
002000 
002010 
002020 
002030 
002040 
002050 
002060 
002070 
002080 
002090 
002100 
002110 
002120 
002130 
002140 
002150 
002160 
002170 
002180 
002190 
002200 
002210 
002220 
0022 JO 
002240 
00223 0 
002260 
002270 
002200 
002290 
002300 
002310 
002320 
002330 
002340 
002J50 
002360 
002370 
002380 
002390 
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>  CONTINUE 

8DN-2 . 726E-04 
SDN- 1. -SDN 
ML  MAX-0. 

WLMIN-100000. 

AMPMAX -0  * 

PO  194  1.^1  »LL 

I- IZ<L)M 

II- UZU  Ml 

I  FORMA  r ( 1X» 1 6  .  F 1 0  • 1 .414.  *  10.3»€15.4.?I4.F10.3»E15.4) 

RR-1 .0-G*R< 1-12 ) 

FK-FC  < I » 1 1 ) *RR tRR 
PO  160J-1.JKAX 
AE-Af S< I. J)/EMI( I»J) 

3  SR*  J-  I 
?  SF-J 

SNVF-SNV(I,II, J, J*1 ) 

SNYR-SNY < I > 1 1 • J» J-l ) 

T  3P-»SNYP*SNYP*SNYP 

r3R-SNYR*SNYR*SNYR 

PL  AF(J) -SDN* 1 OOOOOOOO/SNYP 

PI  AR <  J )  “SDN*  1 OOOOOOOO/SNYR 

MJR< J)-J-2 

MJP< J)-J 

HJQf 

PIP< J)*SP«A£*FK*T3P 
PIR< J>-SR*A£*FK*T3R 
IF<l.N€.  IlK'O  TO  173 
JF(PIP(J). OT . AMPMAX ) AMPNAX-B IP ( J  > 

IF(6JR<J).GT. AMPMAX ) AMPMAX-P IR  <  J  > 

IF(PIF( J).LT. .001*AMPMAX)00  TO  1725 
IF  <  PLAP <  J) ,OT.WLMAX)WLMAX-BLAP< J) 

IF  <  PLAP( J) .LT.ULMIN)WLMIN-PLAP( J) 

5  CONTINUE 

IF ( PIR<  J ) ,LT. . 001  * AMPMAX  > GO  TO  173 
IF ( PL AR ( J  ) . OT . WLMAX ) ULMAX-BL AR ( J) 

IF( PLAR<  J) .  LT. ULMIN )ULMIN"BLAR( J) 

<  CONTINUE 

1F( J. EQ. 20 > WRITE! A. 1901 >  IPATF.  T.  IZ<1_  >  >  IIZ<L>  »MJO< J)>MJP< J) >BLAP< J 
1 . PIP  <  J) »MJQ( J) *MJR<  J) *BLAR( J> »  PIR<  J) 

3  CONTINUE 

IF<  MOO . EO • 1 ) WRI TE  <  7  >(IPATE»T»I7<L)»IIZ(L)  >MJO ( J) >HJP<  J) *  BLAP 

1 J>»BIP< J)»MJO( J) .MJR< J> >BLAR< J) >BIR< J) . J-l , JKAX) 

IF(L. EO. 1. AND. T. EQ. 500. )URITE<6» 1901) (IPATE*TfI2<I.>*IIZ(L)rMJ0(J) 
1 M  JP  <  J ) » Bl  AP  <  J  >  >  B I P  <  J ) »MJO( J) . MJR( J) »PLAR( J) >BIR( J) >  J-l > JKAX) 

I  continue 

IF ( IINU. EG. 1 )CALL  INUEN < l L > JKAX ) 

IFdSPECP.EQ.DCALL  SPECVC ( ULMAX . UL MI N» LLr  JKAX.  JEX.  T ) 

T -T  +  DEL  T 

IF<  T . LE . TMAX  >00  TO  123 
PACKSPACF  7 
PRINT  IVJlf  JKAX 

READ  (  7>  (ID.T.Il.  I2*HJ0<  J)  >M  JP<  J )  »  PI  AP(J)»BIP(J>.MJG<J>» 

IMJR< J).PLAR< J>»bIR< J) > J-l > JKAX) 

WRI TE< 6. t 901 ) (ID.T.Il. 1 2 » MJQ (J) « H  JP<  J) » BLAF'<  J  > »BIP<  J) » M JG( J  >  > 

1HJR <  l>  «Pl  AM  J>.BIR<  J) . J-l .JKAX) 

REWIND  7 

IF ( ICON. EO.l >CALL  HCON< TPCGIN. TMAX . DELT > 

5  CONTINUE 

IF ( ISPECP .EO.l )CALL  PL0TV(4> 

STOP 

END 


SMALL  SUBROUTINES 


FUNCTION  SLOF<A#SO> 

IMP ORFT I CAL  SPECTRA 
SUM  'A*A/(EXP( A/SO> - 1  .()> 

RF  f URN 
FNP 

func r :  >n  snm  t  •  if  > jj> 

:*1M»  NS  I  UN  TC.<2>  *W(2»  <  :*»  »WY<2)  .pt  12)  .ALTAI  ?) 

COMMON  f  •  ,.4f  UX.UY.IU  .AM  ' 

v  r  : 

w  - 1 1  - : 

AJ  .  I 
A  J  «-.».)  1 

!?*»«'  1  >*(UtO.’*  W**  1  »*<‘,»0.‘j)*<V  +  0.5> 

>..*  wi  * » # r vuf  j.*, i  <;»•*<•  -v »  >.5> t<vv*o.5> 

Ml  M  <  1  »  A<  »  AC  |  )  t 1  «/IO.S  * 
f«.  -M  .  2  I  Al  *  A  .'>*<  VV(0.*> 

Dl  T.OiM*  f  J  (I  >*f«i  <  1  »  (UtDItJD) 


002400 
002410 
002420 
002430 
002440 
002450 
002460 
002470 
002460 
002490 
002500 
002510 
002520 
002530 
002540 
002550 
002540 
002570 
002560 
002590 
002600 
002610 
002620 
002630 
002640 
002650 
002660 
002670 
002660 
002690 
002700 
002710 
002720 
002730 
002740 
002750 
002760 
002770 
)  002760 
002790 
002800 
<002610 
002820 
>002630 
002840 
002850 
002660 
002870 
002660 
002690 
002900 
002910 
002920 
002930 
002940 
002950 
002960 
002970 
002960 
002990 
003000 
003010 
003030 
003040 
003050 
003040 
003070 
003060 
003090 
003100 
003110 
003120 
003130 
003140 
003150 
003160 
0031 70 
003  ISO 
003 1^0 
003  •  '0 
00 <710 
003770 
003230 

0  0  <  .’  TO 


D2-4.0«BC<2)tBC<2)«BC<2>/<W<2>tW<2>  > 
AA»AJJ*<AJJM  .0) 

A«AJt(Aj*1.0> 

Fl«Bl*A-PXtAtA 

F2-B26AA-D29AAFAA 

T1*TE<  1  HGiiF  t 

T2«TE<2)+G?4f  2 

SNV-T1-T2 

RETURN 

ENP 


SUBROUTINE  HCONC Ti * T2* T3 ) 

DIMENSION  L<75)*M(75)'A(75>*B<75)fP(7S>*N<73>»C(?5> 

DIMENSION  Y(B01) 

REWIND  3 

DO  100  IFAS-t'l 

REWIND  7 
IV-0 

ST-4800. 

ITI-TI+.001  $  IT2-T24.001  %  IT3-T34.001 

DO  70  IT-irit JT2»XT3 

5  R£AD<7MIP*T*IXrIY*LU>*M(J>*A<J>,B<J>*L<J>*N<  J).C<  J>*D<  J>*  J-l  * 
IF(  UX-IY).NE.IU)0O  TO  3 

DO  6  K-1.001 

6  Y  <K ) *0* 

DO  II  LL*1*4 
1FCLL .£0.1)00  TO  7 

READC  7)<ID'T*IX»lY'L(J>'M(J)rA(J)*P(J)*L(J)'N(J)*C(J)»D<J>*J=l . 

7  IF((IX“IY>  »NE. IV)  GO  TO  110 
ITT»Tf .01 

IF< I TT .NE . IT  >G0  TO  110 

DO  II  1*1 f 2 

DO  11  J»1 *73 

IF(  I . EQ.2)A(J)»C(J) 

1F(1.E0.2>P(J)*D<J) 

Kl-4.t(A< J)-ST)-15. 

K2-K1+30 
IFIK1.LT. 1)K1=1 
IF<K1 . GT .8OO)K1-0OO 
JFCK2.LT » 2)K2*2 
JF(K2.GT .801 )K2=0OJ 
DO  11  K-Kl *K2 
VK*K-l 

DLAM^AC  J>- .254VK-ST 
Y<K)*V(K>4B<J)*SLITFN<IPAS*DLAM> 

11  CONTtNUC 
YH-O. 

DO  12  K-1.801 

IF  <  Y<K ) *GT , YM)KM*K 

JF< Y(K> .GT. YM> YM*Y<K> 

12  CONTINUE 

DO  13  K*1*0O1 

13  YtK)»Y(K>/VM 

WRITE < 3 > IV* T *  I PAS* KM* (Y(K) *K- 1*001 ) 

PRINT  t4*IV*T* I PAS* KM 

14  FORMAT ( IX* 16HSEQUENCE  DECT  U-.I2.7M*  TEMP* *F6 . 0* 1 3M  PEG  K*  IPAS 

1J4.30M  BANDNEAD  AT  STEP  MUH*M) 

70  CONTINUE 
ENDFILE  3 
PRINT  71 

71  FORMAT <2M1*' 

100  CONTINUE 

00  TO  120 

110  PRINT  111 

111  FORMAT (IX* 26MWR0N0  ORDER  ON  DISK*RETURN) 

120  RETURN 

END 


FUNCTION  SLITFNCN'DLAM) 

DIMENSION  FNf 3* 100) • AMUL <3>*LIH(3> 

DATA<FN(1*K)*K>-1* 20 )/5.34*S. 50*3.40' 5. 22* 4.90 *4.56*4. 12*3.64.3.: 
12.76*2.30*1 .06*1.42*1 .00*. 74 *.48 *.20 '.16*. 06*0./ 
DATMFN(?*h)»h-l  *14)/ 7. 93* 7. 83*;  . 37*6 .32*3. 74*4.97*4 . 2* 3. 37* 2. 5 
U.74*l. *.32*. 07*0./ 


003250 
003260 
00J270 
00 32BO 
003290 
003300 
003310 
003 J 20 
003J30 
003340 
003360 
003370 
003300 
003390 
003400 
003410 
003420 
003430 
003440 
003430 
003460 
003470 
003480 
001490 
003500 
003510 
003520 
75)003530 
003540 
003530 
003560 
003570 
00 3 5 80 
7T.  >  0035V  0 
003600 
003610 
003620 
003630 
003640 
003 ASO 
003660 
003670 
003680 
003690 
003700 
003710 
003720 
003730 
003740 
003750 
003760 
003770 
003780 
003790 
003000 
003810 
003820 
003830 
003040 
003850 
003860 
**  003070 
00.3880 
001890 
003900 
00  1710 
003920 
00  iv JO 
0  O  3  V  4  0 
0039J.0 
003960 
003V 70 
003980 
004000 
004010 
004020 
004030 
004040 
004  050 
004060 
22.004070 
004000 
4.  004090 
004 100 
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DATA<FN<3,K),K-l,10>/0.0,8.6,0. I, 7.4,6. 6, 5. 0,5. 0*4. 2.3. 4,2.65.  004110 

11. 0,1. 2,. 75* .4,. 2, .1* .05,0./  004120 

DATA  LIN/20* 14*18/  004130 

DATA  AHUL/. 0794226* .6, .4042/  004140 

C  FOLLOWING  CONSTANTS  FOR  ALADDIN  II  CLARA  I  HOPE  APR  1972  004130 

DATA  LIH/17/  004160 

DATA  AHUL/. 10660/  004170 

PATA<  FN< 1 *K ) t KM • 17) /89 .0*89. , 04 . , 77. , 60 . 4 *39. 4 ,50. 0,  004100 

142. *33.0,25.2,16.4,9.6,5.4,2.0*1 .4# .6*0./  '  004190 

PATA  LIH/13/  004200 

DATA  AHUL / . 1019/  004210 

DATA  (FN< 1 *K>  *K-1 , 131/227. ,224 . ,209. r 181 . , 153. *  123. *90. , 70.  004220 

1 .42. .  16. . 2. .0. »0. /  004230 

C  ALADDIN  74  SLIT  FN  004240 

DATA  LIH/30/  004250 

DATA  AHUL/. 15/  004260 

DATA  <FN< 1»K) *K-1 ,30>/703. * 701 . *609 . ,671 . .644. *616. .307. .550.  004270 

1.329..  500.. 472. .443. .414. .383. .356.. 327.. 298.. 270.. 241.  004200 

2.212. . 184. .135. .126. .97. .68. .41 ..71. .9. .1 . »0./  004290 

C  POST-ALADDIN  JAN  75  SLIT  FN  AL  BURNER  004300 

DATA  AHUL/. 066175. .081125/. LIH/44. 36/  004310 

DATA  <FN( l.K) .K-l .44)/. 0949, .0940, .0945, .093. .091 , .0885, .0059, .083004320 

14..  0809. .0704. .0758. .0733. .0708. .0602. .0657. .0632. .0606. .0501 . .055004330 

26. . 0531 . .0505. . 048. . 0455 .. 0429. .0404. .0379. .0354. .0320. .0303. .0278004340 

3. . 0252. .0227. . 0202. . 01 76. . 0151 .. 01 26. . 01 01 . .0075. .0050. .0027. .0013004350 

4. . 0007.0. .0./  004360 

C  POST-Al ADDIN  JAN  75  SLIT  FN  THA  RELEASE  004370 

DATA  < FN< 2, K)*K-1* 36)/. 1176,. 1173, .1156, .1120, .1003, .1045, . 1000, .0004380 

197. . 0933. .0895. .0050. .082. .0703. .0745. .0708. .0670. .0633. . 0596. .055004390 

28. . 0521 . .0483. .0446. .0400. .0371 . .0333. .0296. .0250. .0221 ,* 0103, . 014004400 

36..  01 00. .0071 . .0034. .0005.0. .0./  004410 

C  AEOLUS  JAN  1975  004420 

DATA  AHUL/. 15345005/, LIH/29/  004430 

DATA  (FN(1,K) , K-l, 29) /142. 5, 140. ,134. ,127.8,121 .7*115.6,109.4,103.004440 

12.97.2.91 . . 84.8.76.6.72.5.66.3.60.2.54.2.48.1.42.0.35.8.29.4.23.2.004450 

217.0,11 .5, 6.5, 4. ,2. ,.7,0. ,0./  004460 

DS"AB5<  PLAN ) /AHUL (N) 41 •  004470 

IDS-DS  004480 

IF(IPS.0E.LIH(N)-1 >00  TO  20  004490 

CDS'- 1  PS  004500 

RDS-DS-DDS  004510 

SL I  TFN=*  ( 1 .  -RDS  > *FN<  N,  IDS )  f  RDS4FN  (N,  IDS  +  1  >  004520 

00  TO  30  004530 

20  SLITFN-O.  004540 

30  RETURN  004550 


END 


004560 


50 


PROGRAM  MArCH'OUTPUT. TAPES* TAPE6. TAPE2. TAFf J)  008?  10 

DIMENSION  N8F  < • >  » NSTR  <  0 ' •  MICRON  <  8 ) . FPL  I  ST  <  J  >  .AMUL  <01  008??0 

DIMENSION  AMBDAP I  3  *  008?  JO 

DIMENSION  WAC <  20  > . WA* ( 20 )  008240 

COMMON  V ( 300  >  *  X<  300  >  *  YW<  300  > .  Z  < 80 1 >  * /W<  801 >»TI19>.  DUN <  801 1  006230 

INTEGER  ON  00U260 

DATA  A9KY/10MBLUC  S*  1  /  008?7O 

C  DATA  AMUL/.  3971 13.  .397  1  1.1.  .J2408.’,  .324082.  .324002.  .340103*  0082BP 

C  1  1  .62444*  1 . 62444/  008/90 

DATA  AMUL/. 2212. .27|2. .271?. .271?. .2/12. .2712. .?/|2,  .2712/  008300 

DATA  N3F/1 »1»1'1» 1*1.1. 1/  008310 

C  DATA  NSTR/60.40. 183.278.67.231 .291/  008320 

C  OATA  MStR/224. 61 , 195, 30.203. 112.262. 14  \/  006330 

C  ALADDIN  II  ClARA  AND  HOPE  008340 

C  DATA  NSTR/ 170.1 49* 144*? 30*0. 0*0.0*/  008330 

C  GOLIATH  NOV  71  008340 

DATA  AMUL/. 4074. .4078*0. *0. *0. *0. *0. *0./  008370 

DATA  NSTR/46. 1 1 7 . 0.0. 0.0. 0 .0/  008380 

DATA  HICR0N/1J5. 1 J5. 1 15* 1 7A»1 76. 17A.I 74.174/  008390 

C  ALADDIN  74  008400 

DATA  NSTR/79  *  79 . 79*  79. 79*  79 . 79*  79/  008410 

DATA  AMUL/. 3035' .3035. . 3053. .3055.0. *0. .0. »0./  008420 

C  POST-ALADDIN  JAN  75  0084.30 

DATA  NSF/1 » 1 .2.2* ?.?»?» 2/  008440 

DATA  NSTR/204 ♦ l *  224  »1*1*1*1»1/  008450 

DATA  AMUL/. 2647.1... 3245. 1 .  »  1 .  .  1 . » 1 . .  1 » /  008460 

DATA  AMBDAS/4824 . 3/ * AMBDAF /4886 » 3/ » AHBDAP/3B4042 .5/  008470 

DATA  AMUL/. 3069* .3069. .3069. .3069* .3069* .3069. .3069. .3069/  000480 

DATA  NSF/1 *1. 1  *  1 . 1 . 1 • 1  *  1/  000490 

DATA  N6TR/290. 143. 244. M3. 192. 143.1.1/  008500 

C  USMR  OCT  77  008510 

DATA  NSTR/8* 180/. AMUL /0*. JO 10/  008520 

REWIND  2  008530 

REWIND  3  008540 

REWIND  5  000550 

C  IF  NUMBA.GE .6  A  CRT  PLOT  IS  MADE*  OTHERWISE  A  PEN  PLOT  008560 

C  IF  MUMBA*0  NO  PLOT  IS  HADE.  008570 

CC  IAP  IS  PARAMETER  FOR  SCALE  ON  APLOT.  000580 

C  IF  SN»0.  NO  SMOOTHING  IS  DONE.  008590 

C  SM  IS  SMOOTHING  PARABOLA  HALF  WIDTH  IN  ANGSTROMS.  000600 

C  IF  STPL  -  0.  NO  STRIP  PLOT  IS  MADE.  008610 

C  STPL*1.  MAKES  STRIP  PLOT  ON  OUTPUT  WITH  APLOT.  008620 

C  IF  STPL-2.  THE  STRIP  PLOT  IS  MADE  BUT  NO  OTHER  PROCESSING  IS  DONE.  008630 
C  IF  STPL*3.  ONLY  STRIP  PLOT  IS  MADE  AND  VALUE  READ  IN  AS  PFAK  IS  USED00B640 
C  FOR  NSTR.  008450 

C  IF  STPL-4.  ONLY  STRIP  PLOT  IS  MADE.  PEAK  IS  USED  FOR  NSTF  008660 

C  AND  SHOT  IS  USED  FOR  SMOOTHING.  008670 

C  IF  OW*QUICK  PLOT.  QUICK  CRT  PLOT  OF  FILES  N1  *  N2.  OF  DATS-  TAPE  IS  008680 

C  MADE  BUT  NO  OTHER  PROCESSING  IS  DONE.  000690 

READ ( 5 . 2 1 0 ) OW . N 1 . N2  *  S  T  PL . I AP , SM . NUMB A l . NUMB A2  008  700 

210  FORMAT<4X»II.6XrI2.3X.I2.7X.F2.0»3X*I2.5X,F2.0*6X*Il»5  .11)  000710 

PRINT  211.0W.N1 . N2.STPL* IAP » SM . NUMBA 1 .NUMBA2  008720 

211  FORMAT  < IX*  3H0W*  *  1 1  *  1 X . 5HF ILFS  * I2.1Xf?HT0*I2»7M*  STPI *.  2.0*  008730 

11X»4HIAP».  I2.5H.  SM=*F?.0*6H*  PEN=* *  T  l .  IX . 4MCRT- *  1 11  008740 

NUMBA»NUMBAM7*NUMBA2  008750 

IF ( OW .NE .1)00  TO  1  008760 

CALL  AL2<N1 .N2.NSTR )  006770 

GO  TO  100  000700 

1  CONTINUE  000790 

PRINT  201  008000 

201  FORMAT <  02H1FILE  NO.  SCAN  NO.  SCANS  SUMMED  TEMPERATURE ( DEG .  K) 0080 10 

1  STD  DEVIATION  DESCRIPTION. 9X . 9MPEAK  SHOT)  008820 

XPEAK-O.  0088.30 

DO  3  K*i .801  000P4C 

3  ZW(K )*4800« ♦ . 25BFL0AT  <K-l  >  008850 

DO  5  L«lfl9  008860 

5  T (L ) “504LF  tOO  OOHH/O 

7  READ  (5»20: ) NF .NSC.N3UM. ALPH1 . ALPH2.PE AK . ALPH3  008080 

IFCEOF (5) ) 100.B  000090 

202  FORMAT  <9X*T?.10X*I4.5X»I3»9X.A10.A9rFA.?»A5>  008960 

0  IFtSTPL.LT.l. 5 > CALL  ZSTORE<-NSF(NF) )  0089(0 

CALI  I SNN  <  NF • M  .  I DUMMY )  009920 

000930 

NSUM=MAXO<  NSUM. 1 )  000940 

DO  10  N* l*30C  000950 

10  V<N)*1 .  008960 

N9TRNF -NSTRfNF  000970 

C  AEOLUS  MARCH  1975  006980 

TF<NSC.GT.58>NSTRNF»NSTR(NF+2»  008990 

IF ( NSC . G  r . 69 )NSTRNF=NSTR( NFE4 )  009000 

IF ( SI TL  »6T . 2. 9 )NSTRNr  PEAK  009010 

DO  2?  N- t .NSUM  009020 

N4  *.3006  ( Nf  NSC  >  ENSTfcNF  -601  009030 

DO  20  K  *  1 . 300  0090  40 
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CAU.  ISNN(R)N4*-1 # IWNJ) 

20  Y(K)-Y(K>  fTKAYt IUN3) 

2?  CONTINUE 

IT (STPL.tE. .9)00  TO  25 

ANSUM-NSUM 

DO  23  K- 1  * 300 

23  V(K> -Y<K >/AN5UM 
PRINT  201 

UR1 TE  <  6*  201 ) 

PRINT  204 #NF  tNSC  t NSUM* AL PHI » ALPH2 *  PEAK  r ALPH3 
UR  I rc  <  6 * 204 ) NF  I NSC  *  NSUH . AL  PHI  »  AI.PH? • PEAK  *  ALPH3 
IFISTPL .GE *3*9)00  TO  235 
CALL  APLOT ( 300  # I AP ) 

IF (STPL  *GE • I *V)GO  TO  7 
235  00  24  K« 1.300 

24  Y<KMY<K>*ANSUM 
75  NSFNF ^NSF<NF> 

C 

C 

C 

C  FUR  INTERACTIVE  ADJUSTMENT 
INPASS  *  1  t  GO  TO  27 

c 

INPAS5=0 

76  NUMBA-0  4  INPASS- INPASS41 
PRINT  1002 

267  READ  1001 » APUT 

IF ( APUT • NE . 0 • )G0  TO  264  %  NUMBA-10  f  00  TO  262 
264  PEAK=APUT 
PRINT  1003 
REAP  IOOIfAPUT 
IF ( APUT . NE . 0 . ) AMUL ( NF ) “APUT 

1001  FORMAT  <F13. 2) 

1002  FORMAT <4  PEAK-*) 

1003  FORMAT (*  AMUL“*> 

27  CONTINUE 

C 

C 

c 

c 

00  30  K»1 f 300 

YU(  K  )  “AMBDAP <NSFNF  )  -f  AMUL  <  NF ) * < PEAK-FLOAT (K)  ) 

IF< INPASS. EG. 1 ) Y(K ) ~FRES<  Y(K ) f YW<K ) ) 

30  CONTINUE 

IF  I TS“PEAK- <  AMBDAF -AMBDAP <  NSFNF ) ) /AMUL ( NF ) 

I r I TF  =PE AK 4 <  AMBDAP ( NSFNF  > - AMBDAS ) /AMUL ( NF ) 

1FITS-MIN0T IF  ITS » 275 ) 

IFITF“MAXO( IFITFf 1 ) 

IFIT5-MAX0(IFITSf1 ) 

IFITF-MIN0<IFITF*273) 

IF(STPL.GE.3.9)IFITS-1 
IF (STPL . GE  . 3.9) IF  I TF-300 
IF<STPL.LE.3.9>G0  TO  37 
CALL  API.  OT  (  300*  I  AP ) 

IF< STPL . GE . 1 . 9 )G0  TO  7 
37  IF< ASKY . NF . ALPM1 )G0  TO  50 
DO  40  K - 1  *  300 
40  X(K)=Y<K) 

XPEAK=FEAK 

PRINT  204  * NF * NSC  * NSUM* ALFU1 » ALPH2* PEAK* ALPH3 
204  FORMAT  <3X* I2*OXfI3f 1  OX  *  1 2*  43X  *  A1 Of  A9  *F6 . 1  *  1 X  r  A5 ) 

GO  TO  7 
50  CONTINUE 

CALL  FITCROl  *19*  IF!  TSfIFUFfTF*  TE  fFM  *  WAC  *  UAK  *  0 .  •  XPEAK  *  AMUL  <  NT  )  ) 
IDLIST(1)=NF*IDLIST(2>“NSC*IDLIST(3)-NSUM 
PR I N  T  203  *  NF  *  NSC  *  NSUM  *  TF  *  TE  *  ALPM 1 r  AL  PH2  *  PE  AK » AL  PM3 
703  F0RMAT<3X»I?.flX*I3. IOX* 17* I4X*F6. 1 »I3X.F5. 1 f5X»AI0fA9fF6. 1 • 1 X»A5) 

C 

C 

C 

C  FOR  INTERACTIVE  ADJUSTMENT 
C  T F ( NUMB A . NE . 1 0  > GO  TO  76 

t: 

C 

c 

IF (NUMPA.EG.O)GO  Tn  7 
CALL  SM0T<YII*Ytl.773fSM) 

CALL  NEUPI  TdDLIST  *1  *3.001  *  l».*/5*TF  ♦  T  F  .MI  CRON  <  NF  )  *  1 9  *  UAC  » UAK  * 

1 0  *  NUMB A  *  SM  *  XPf  AK  *  AMUL ( Nl  )> 

GO  TO  / 
too  roNriNtM 

IF  (NUMB  A  ,  I  (J.  |  .OF  .NIIMBA.'  U.  7>CAi  I  ENPH  T 
Of  01 


009050 
009060 
009070 
009000 
009090 
009100 
009110 
009170 
009130 
009140 
009150 
009160 
009170 
009180 
009190 
009200 
009210 
009220 
009230 
009240 
009250 
009260 
009270 
009280 
009290 
009300 
009310 
009320 
009330 
009340 
009350 
009360 
009370 
009300 
009390 
009400 
009410 
009420 
009430 
009440 
009450 
009460 
009470 
009480 
009490 
009500 
009510 
009520 
009530 
009540 
009550 
009560 
009570 
009580 
009590 
009600 
009610 
009620 
009630 
009640 
009650 
009660 
009670 
009680 
009690 
009700 
009710 
009770 
009730 
009740 
009750 
009760 
009770 
009700 
009790 
009000 
0098 10 
009020 
009030 
O09840 
009030 
009060 
OO°0  70 


f  NH 


009HP0 
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on  non  n  n  nnnnn  n  nnnnn  nnnnr 


12  FRF9 

009900 

009910 

009920 

FUNCTION  FRES(F.FLAMD> 

009930 

009940 

009950 

FL-FL AMP-4800. 

009960 

USED  WITH  HAWAII  DATA 

009970 

FRES-F9<  .  930753* .  001 19154FL- .  00000504*»LtFL  > 

009980 

USED  WITH  ALADDIN  AND  GENIE  DATA 

009990 

oioooo 

010010 

FR£S*F*< 1.4. 00049034FL ) 

010020 

RETURN 

010030 

010040 

END 

010050 

13  FKAY 

010070 

010090 

010090 

FUNCTION  FKAYIKAYPRM) 

010100 

010110 

010120 

ALADDIN  74 

010130 

DATA  A/98iB.9/,B/.0004572/.T/.0011418/ 

010140 

POST-ALADDIN  JAN  75 

010150 

DATA  A/1 ./.B/1.OO780/.T/. 0011621/ 

0)0160 

WSMR  OCT  77 

010170 

010100 

010190 

FKAY«KAYPRM 

010200 

RETURN 

010210 

FK-KAVPRM 

010220 

010230 

010240 

USED  UP  TO  APR  73  ON  ALADDIN  11  DATA 

010250 

FKAY-FK/U.-FK/950*  > 

010260 

SLIGHTLY  DIFFERENT  NO.  FOR  GOLIATH  NOV  71 

010270 

FKAY-FK/ ( 1 ♦ -FN/949 .217) 

010280 

AEOLUS  MARCH  L975 

010790 

DATA  A/2699 . 0/ . 8/ . 0002055/ . T/ . 00 1 0493/ 

010300 

FKAY*FK/< 1 . -FK*T)-B*A 

010310 

RETURN 

010320 

END 

010330 

c 

c 

C  14 

C 

C 


FIT  POISSON 


SUBROU TINE  FIT<N»NT.N1 »N2»TF»  TE  »FM. « ACC  *  U AKK t S H , XPE  AK  * AMUL ) 


010350 

0103A0 

010370 

010300 

010390 

010400 

010410 

“010420 


WAVELENGTH  ARRAY  SYNTHETIC  SPECTRA  N  El  EMENTS 
AMPLITUDE  ARRAY  NT  BY  N 


ZU 

1 
H 

T  TEMPERATURE  ARRAY  NT  ELEMENTS 
u  ▼  miict  nr  i  rqc  tuam 

YU  MEASURED  SPECTRA  WAVELENGTH  ARRAY  NY  ELEMENTS 
Y  MEASURED  SPECTRA  AMPLITUDE  ARRAY  NY  ELEMENTS 
NY 

TF  BEST  FIT  TEMPERATURE 
TE  ERROR  BARS  ON  TEMPERATURE 
X  -  BACKGROUND  ARRAY 

NOTE  THE  SYNTHETIC  ARRAY  MUST  BE  IN  CONSTANT  INCREMENTS 


010430 
010440 
010450 
010460 
010470 
010400 
010490 
010500 
010510 
010520 
010530 
010540 
010550 
IN  WAVELENGTO 1 0560 
010570 
010580 

- 010590 

010600 

DIMENSION  UACC(20>  »UAKK ( 20 ) » AA<  20 ) » AMI  20 ) » F( 19)  010610 

COMMON  Y ( 300  >  * X < 300 ) • YW(300) *  Z ( 00 1 >  * ZU< 001 > »  T( 19 > » DUH (001)  01 0620 

AN=N-l  010630 

I =N1  010640 

CONTINUE  010650 

AK-AN4<YU<I)-ZW(1))/(ZU(N)“ZW(1>>41.  010660 

IF(AK»GE< 1 » .AND.AK.LT * ANfl . >00  TO  20  010670 

IFCAK-l. >13*13.14  010600 

13  N1*I41  010690 

GO  TO  20  010700 

14  N2-I-1  010710 

GO  TO  21  010720 

IF <  I  »Gf . N2 >00  TO  21  010730 

I » 1+ 1  010740 

GO  TO  10  010750 

21  D~N2-N1+1  010760 

YMAX=0.  •  YMIN-1.E+O0  t  BYC*=0.  010770 

DO  40  I *N1 *N2  010780 

YMAXaAMAXl ( YMAX* Y ( I ) )  010790 

YMIN^AMINl <  YMINr  Y< I ) )  010000 

IF<Y< I ) .  L  T  •  5  • )  GO  TO  401  010010 

VI 1  *-Y<  I  > 41  •  010020 

PVr-HYCf (YI1-.5)*A|  OG<VI| >-YI141 ./(12.4YI1 > -1 . /< J60 . *Y 1 1**3 )  010830 

Hl./< 1260.«Yllt*5>- 1 ./< 1600. *YI1 **7>  010040 

GO  111  405  010050 

1  K)Y  I  010860 

KY  >1)4.5  010070 

IT (KY.I  0. 2>M>  V»2  010000 

If  Oi.dl,  J'M  TIOA'-O 

I)  fM  .(  U.41M*.  i  -24  A 1 0900 

IM*  r.rr>.5)fcfcYeI?0  010910 

YM\Y -f  l  OAT  <KKY )  010920 

BYC - BYC4AL0G  <  YKK Y  >  010930 

5  Y  IP- AMAX1  <«01»Y(I)>  010940 

DYC-RVC- .5*AL0G< Y I*  ;  010950 

CCNTINUe  010960 

U4K* YMAX-YMIN  •  UAC»YMIN  010970 

UO  5"  2-1 »NT  010900 

CAM  7S  f  ORE  (  J  >  .'10VVQ 

LUOf^O  7M0OO 

A  7-0.  i  AY/r 1  0.  •  AYfl. -0.  *  AY72D7  O.  ♦  AV7n/=0.  4  AYP7*0.  OUO10 

incip-innrn  oho/o 

I  r  U  UOf  .  G I  •  50  )  T  (  J  > «  1  •  I  )  1  1  ••'MOT) 

U  <  0U»  .G  1.50)00  10  5  »  011OI  - 

1*0  46  I  *N1  *N?  * 1 1 01  •' 

AK -AN»< YW< I >  ZU< I) >/</W  N>  7U.1))41.  • 

R' AK  0110/0 


cr  -k 

O 1  .000 

*>K  »A>  l> 

oi  iovo 

i /  -  1 i .  AR  •  /<* ) »/<K4f  >*AK 

01 J lCrt 

A/  A7*///Il 

•>11  1  JA 

0A  t  UAK4/7  fUAT 

01 1 1/0 

b»n  I  >/  (P*BA  > 

111  lo 

H}‘  -0/8.> 

1 1 ;  ao 

AYU-AYD t B 

01  1150 

AY PC*  -A  * |«? 1  BH 

o  !  1 160 

i*  !*•/ ; 

01  ’  1  'O 

At/p  a-.iup 

••»  1  1«0 
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•••/BA  011190 

AYZD2-AYZD24*  011200 

B-MZZ  011210 

AYZ2D2-AY72D2fB  011220 

40  CONTINUE  011230 

DET-AYZ2D24AY02-AYZD2442  011740 

FNOT«AZ-AYZD  4  GNOT-l.-AYO  011250 

PELK- <  ONOT  4AVZP2 -FN0T4AYP2 ) /OET  0112*0 

D£LC*<FN0T*AYZD2-0N01»AYZ2B2>/D£T  011270 

WAK-UAK+PELK  4  UAC-UAC1 PCLC  011780 

IF<PELK.GT.1.E-10«UAK.O«.BELC.OT.!.E-10*UAC>GO  TO  44  011290 

MACC(J>-NAC  011300 

UAKK(J)-MAK  011310 

IF1UAC.LT.0. .OR.WAK.LT.O. >F< J>«1 .E+13  011320 

IF 1 WAC  .1  T  .0. • OR. WAK . LT . 0 » 100  TO  SS  011330 

F<J>*0.  011340 

00  50  I-N1.N7  011350 

AK«AN*(YU<I)-ZU<1>  >/(ZW<N)-Z«<  11)41.  01  1 360 

K-AK  011370 

CK-K  011380 

AK-AK-CK  011390 

ZZ"U.-AK)4Z<K)4AK4Z(K41>  011400 

XI-1.  011410 

IFIXPEAK.LT. .001 )G0  TO  48  011420 

FK-XPEAK  - 1  YU  <  I  )  -484 1 . 34 1  / AHUL  011430 

H»FK  011440 

ACO-H  011450 

ACG-FK-ACG  0114*0 

X1*ACG»X<NE1)4<1.-ACG>*X<H)  011470 

48  B-UAK  4ZZ EUAC  011480 

B-B~Y(I>*AL0G<8>  011490 

SO  F<J)«F<J>+B  011500 

55  F ( J) “F  <  Jl+BYC  011510 

201  FORNATl 1X.SE11 .4)  011570 

FH-10.EE10  01 15J0 

DO  *0  3-1 .NT  011540 

IF1F1.J)  .LT.FN)  JT-J  011550 

IF(F(J).LT.FN1FN-F<J)  0115*0 

40  CONTINUE  011570 

IFlJT.EO.l.OR.JT.EO.NT)  GO  TO  71  011580 

Xl-TCJT-1)  011590 

X2-T  <  JT )  011*00 

X3=T(JT+1>  011*10 

PO  70  NN-1,3  011*20 

GO  TO  <*1.*2.*3).NN  011*30 

*1  YJ-F1JT-1)  011*40 

Y2*F< JT)  011*50 

Y3-F1 JTE1 )  011**0 

GO  TO  44  011*70 

*2  Y1-MACC1JT-1)  011*90 

Y2-MACC1 JT  )  011*90 

Y3-WACC1 JTtl)  011700 

GO  TO  *4  011710 

*3  Yl—UAKK ( JT-1 )  011720 

Y2-WAKK1 JT  )  011730 

Y3-WAKKC JT+1 )  011740 

44  Y4»(Y2-Y1)/(X2-X1)  011750 

Y5«<Y3-Y2)/<X3-X2)  0117*0 

X4-1X24X2-X14X1 1/1X2-X1 )  011770 

X5-«X3*X3-X2*X2)/(X3-X2)  011780 

AQ*(Y5-Y4)/<X5-X4)  011790 

B0-Y4-A8BX4  011800 

C0-Y1-A0«X1*X1-BG*X1  011810 

GO  TO  <*3»«A.47>.NN  011870 

*5  TF— BO/12.BAO)  011830 

FN-AOBTFBTF ♦B04TF+C0  011840 

RH-FH*<I.E1./D>  011850 

TE-BG4B0-4  *  4AQ4 <  CO-RH )  0118*0 

TE-SQRT ( TE 1 7 ( 2 . * AO )  011870 

GO  TO  70  011880 

A*  IIAC-AO*TF*TFEBO*TFECO  011890 

GO  TO  70  011900 

*7  MAK-A04TF4TF+B04TF4CG  011910 

70  CONTINUE  011920 

00  TO  72  011930 

71  TF-T(JT)  011940 

TE-999.0  011950 

72  RETURN  0119*0 

END  011970 

01 1980 
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c 

012000 

c 

012010 

C  13 

ZSTORE 

012020 

C 

012030 

c 

012040 

SUBROUTINE  ZSTORE(L) 

012030 

COMMON  Y<300),XC300>*YW<300>»Z<801>»ZU<801 > » T < 19 > , DUM10O1 > 

012060 

IF(L.LT.l)  00  TO  16 

012070 

1FU-LL  >91 ,100.6 

012080 

6 

LM‘LL*1 

012090 

00  B  M*LM»L 

012100 

fl 

READ< 3 > IV, TT»MM»KK #  <  Z< J J> » JJ-1 ,801 > 

012110 

LL»L 

012120 

00  TO  100 

012130 

16 

N5- 132670981 22369 

012140 

IF<N4.NE.N5>N2-M3 

012150 

N4-*N5 

012160 

NF  =  -L 

012170 

IF1NF.GT .N2)00  TO  22 

012180 

20 

REWIND  3 

012190 

N2»l 

012200 

it 

N3*N2+1 

012210 

IF(NE.EO.l)  GO  TO  20 

012320 

73 

DO  23  N1»N3,NF 

012230 

24 

READ< 3) IV. TT rMM.KK , (  Z< JJ> , JJ«1 ,001 > 

012240 

IF(E0F<3>>25,24 

012250 

73 

CONTINUE 

012260 

28 

LL*0 

012270 

90 

N2«NF 

012280 

ira-i  >100,6,6 

012290 

91 

BACKSPACE  3 

012300 

92 

LL^l L~1 

012310 

BACKSPACE  3 

012320 

IF<L  .LT.LDOO  TO  92 

012330 

REAP< 3> IV, TT»NM»KK, <  Z< JJ> , JJ-1 ,801) 

012340 

100 

CONTINUE 

012350 

RETURN 

012360 

FND 

012370 

C 

012390 

r 

012400 

r  ia 

SHOT  RECTANGULAR 

012410 

r: 

012420 

c 

012430 

SUBROUTINE  SM0T1X, Y,N1 ,N2,A> 

012440 

COMMON  YD(300> , XD<  300) , YWD(300> , ZDC80I > ,ZWD<801 >  »TD< 19) ,Z<0Ol > 

012450 

DIMENSION  X< 1 ) , Y < 1 ) 

012460 

1FIA.LE. .001 >00  TO  100 

012470 

AS*<  X(N1 >-X<N2> > / (N1 -N2  > 

012400 

ASSADS < AS > 

017490 

II**  56 A/ ASF « 5001 

012500 

NN-24 II-l 

012510 

I I* I I~1 

012520 

NN1 =N1 ♦ I I  t  NN2-N2-II 

012530 

DO  30  I«=NN1,NN2 

012540 

Z(  I  >  <<  I > 

012550 

DO  30  J«1,II 

012560 

1.1*1- J  t  JMIJ 

012570 

>F(Y(IJ>F» < JI> > 

012580 

50 

CONTINUE 

012590 

Dfl  40  I-NN1,NN2 

012600 

40 

V< I >*Z( I >/FLOAT  <NN) 

012610 

100 

fcF.llIRN 

012670 

END 

012630 

C 

012650 

r 

012460 

r  i  / 

I  !»NN 

012670 

r. 

012600 

f. 

01 2ft 90 

SUBROUT INF  I  SNN(N,K , ISN> 

0 l 7 700 

D 1  Ml'  NS  I  ON  T I  (  30  > 

01 2710 

ir<K>20,20,10 

012720 

10 

N6~ 1 3207690! 239 

012730 

IF(N4,NF.N6>  N2*NA 

012740 

N4  N6 

012750 

NF  =  N 

012760 

GO  TO  100 

0127-0 

70 

IF  *  NF -N2 >21 , 11,22 

012780 

?: 

REWIND  2 

012790 

N5»-2V 

012000 

if  <NF  *m.  1  >00  TO  32 

012010 

N.»*  1 

0 ! ; 020 

Nr.-  79 

o  1 4*n  <o 

N5  N.'M 

012040 

M»  30  Nf  '  N3*N» 

01  3ll?.0 

56 


25 

READ <  2*209)  < 11 <NN)  *NNM  *  J0> 

012860 

IF (EOF  < 2) >30*  25 

012070 

30 

CONTINUE 

012000 

00  TO  32 

012090 

31 

IF(M.LT.M5>G0  TO  50 

012900 

IF(N.GE. N5f30>00  TO  32 

012910 

GO  TQ  01 

012920 

3? 

N2-NF 

012930 

N5aN5+30 

01 7940 

DO  4*  N1-N5*N*30 

012950 

N1 1*N1 

012960 

40 

READ( 2  *  209  > ( 1 1 (NN>  *NN»1 *  30) 

012970 

209 

FORMAT  <  3013 ) 

012900 

N3-NU 

012990 

00  TO  01 

013000 

50 

BACKSPACE  2 

013010 

Si 

N3«NS-30 

013020 

BACKSPACE  2 

0t3030 

IF(N.LT.N5>  00  TO  51 

013040 

READ  ( 2 » 209 HII(NN)* NN» 1*30) 

013050 

01 

NN'N-N3+1 

013060 

ISN*II<NN> 

013070 

100 

RETURN 

013000 

END 

013090 
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C 

SUBROUTINE  NEUPLT  ( IDLI8T  *K 1  rK2*N8YP* N1  *N2*  rr  *  TE*  MICRON  *  NT  *  ACC* 
1AKK  * IV* RUMBA *  8M * XPEAK • AHUL ) 

C 

C  HAKES  FANCY  PLOT  OP  SYNTHETIC  SPECTRA  COMBINED  WITH  MEASURED  SPECTRA 
C  I DL 1ST <k 1 >  TO  IDLIST<K2>  IS  A  LIST  OP  SCAN  NUMBERS  TO  BE  PLOTTED 

C  Z(RTfNSYP>  IS  AN  ARRAY  OP  SYNTHETIC  SPECTRA 
C  ZU  IS  THE  ARRAY  OP  WAVELENGTHS  POR  THE  ELEMENTS  IN  Z 
C  Y  IS  AN  ARRAY  OP  MEA8URED  SPECTRA  WITH  NYEl.  ELEMENTS 
C  YU  IS  THE  ARRAY  OP  WAVELFNGTHS  POR  THE  El  EMENTS  IN  Y 

C  rr  IS  THE  BEST  TIT  TEMPERATURE  POR  THAT  SUM  OP  SCANS  BY  UP  IT 

C  TC  IS  THE  ERROR  BAR  ON  TP 
C  MICRON  IS  THE  SLIT  WIDTHE  CODE  IN  MICRONS 

C  T  IS  THE  ARRAY  OP  TEMPERATURES  USED  IN  THE  SYNTHETIC  CURVES 
C  NT  IS  THE  NUMBER  OF  TEMPERATURES 

C  ACC  AND  ARK  ARE  THE  BEST  ADDITIVE  AND  MULTIPLICATIVE  CONSTANTS 
C  IV  SFOUENCE  OF  ALO  USED 
C 

DIMENSION  IDLXSTU) 

DIMENSION  ACC <  20 > *  AKK <  20 ) 

COMMON  Y ( 300 >  *  X  < 300 >  *  YU <  300  >  *  Z ( 00 1 >  *  ZU < 00 1 >  *  T<19 >  *  YP < 801 > 

INTEGER  H 
C  DRAW  AXIS 

IF(NUMBA.EQ.10>00  TO  S 
CALL  SUPAPL<IV*2*1 *NUMBA*0*0> 

GO  TO  6 

5  TID-IDHST(3> 

CALL  PPLOT (1*ZW*Z*IDLIST(1)*IDLI9T<2)*TID> 

C  FIND  NUMBER  OP  BEST  PIT  TEMPERATURE 

6  TT*T<NT ) 

DO  1010  N-1'NT 
DT«ABS(T  (N)-TP) 

IF<DT.LT.TT)LL-M 
1010  IF<DT.LT.TT)TT-DT 
YMAX-O. 

DO  8  M-N1*N2 

0  YMAX-AMAX1 ( YMAXpY(M) ) 

C  IP  NUMBA-7  A  CRT  PLOT  WILL  BE  MADE. 

C  IF  NUMBA-10  A  PRINT  PLOT  UII  L  BE  MADE. 

C  IF  NUMBA-1  A  PEN  PLOT  UILl  BE  MADE. 

C  FOR  PEN  PLOT  PP  NAMED  PENPLOIS  MUST  BE  LOADED. 

C  FOR  CRT  PLOT  PP  NAMED  CRTPLOTS  MUST  BE  LOADED. 

C  FOR  PRINT  PLOT  PP  NAMED  KI fMAPAB  HU8T  BE  LOAPED. 

NDI-NYEL-3 

C  PLOT  SYNTHETIC  SPECTRA 
DO  32  I-WNSYP 
I2-I-1 

IP(ZW<n.OT.YW<Nl)  )0D  TO  33 

32  CONTINUE 

33  DO  34  I-ltNSYP 
11*1 

I F ( ZW  < I ) . OE  *  YU ( N2 ) ) GO  TO  35 

34  CONTINUE 

35  DO  1100  MM-1*5 
L-LLFHM-3 

IF(L.LT.l .OR.L.OT. 19)00  TO  ItOO 
CALL  ZSTORC<L> 

C  SHOT  CHANGES  PRESENT  Z  ARRAY 
CALL  SM0T<ZW*Z*I1*I2*SM> 

IF ( MM . NE . 1 )G0  TO  37 
37  DO  30  IM.N5YP 

Z<I)»<Z<I)*AKK<L)4ACC<L))/YMAX 
Z<I )*AMAX1 <Z<t) * .01 > 

30  ZU)-4.«AL0G10(Z(I)  )>8. 

IF (NUMBA- 10)40* 39*40 

39  CALL  PPLOT <MMM  *  ZU*  Z* 1 *NSYP*T<L ) ) 

GO  TO  1100 

40  IP-*3 

DO  1090  M*  1  *N!»VP 
ZZ-Z(M) 

X3'(ZU  (M)  -  4000.  F200. 4PI.0AT <  IV)  >4.08 
ir(ZZ.OT.0. .OR.ZZ.LT, 0. .0R.X3.LT. 0. )G0  TO  1090 
IP(X3.0r.l2. >00  TO  1 096 
CALL  PLOT (X3*ZZ* IP> 

1090  IP-2 
1096  AL-MM 

C  LADLE  SYNTHETIC  SPECTRA  POR  TEMPERATURE 
CALL  NUMBER! 12.5* Al *.2*T  (L>»0.»-1) 

CALL  SYMBOL  < 13«3*AL* . 2*5HDEG  K*0.*5) 

1100  CONTINUE 
C  PLOf  MEASURE  D  F»l  I.C TRA 
IP  *3 

DO  I0/.9  M-Nl.N? 


013140 
013150 
013160 
013170 
013180 
013190 
013200 
013210 
013220 
013730 
013740 
013250 
01 32AO 
013270 
013280 
013290 
013300 
013310 
013320 
013330 
013340 
013350 
013360 
013370 
013380 
013390 
013400 
013410 
013420 
013430 
013440 
013450 
013460 
013470 
013480 
013490 
013500 
013510 
013520 
013530 
013540 
013550 
0 1 3560 
Ot  3570 
013500 
013590 
013600 
013610 
013620 
013630 
013640 
013650 
013660 
013670 
013600 
013690 
013700 
013710 
013720 
013730 
013740 
013750 
013760 
013770 
0  1 1 700 
013790 
013000 
01 JH10 
01 3H70 
013030 
0131140 
0 1 3050 
013060 
013070 
013880 
013890 
013900 
013910 
013920 
013930 
01.1940 
0 1 3950 
013760 
013970 
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X2«!YU  <H)-4800.4200.»FLOAT<IV>>*.08  013980 

XN-1 .  0139*0 

IFIXPEAK.LT.. 001)00  TO  1035  014000 

FK»XPEAK-!YW  IM1-48A1 .341/AMUL  014010 

I -EX  014020 

ACO-1  014030 

ACO-FK-ACO  014040 

XN-AC0*X<I41>4<1.-ACG>«X<I)  01 4050 

1035  YY"100.*Y!N)/YNAX  0140*0 

yy»amaxi<yy.i.>  014070 

Y3"4 . 4AL0G10  <  YY >  0140B0 

YP!N)-Y3  014090 

1F<X2.LT.0. .0R.Y3.LT,0..OR.Y3.0T.9,)00  TO  10*9  014100 

JF<X2.GT.12. >00  TO  1070  014110 

IF!NUMBA.EO.  10)00  TO  10*9  014120 

CALL  PL0T!X2.Y3»IP>  014130 

1049  014140 

1070  IFINUMBA.NE. 10)00  TO  1071  014150 

CALL  PPL0T!7.YW.YP.N1.N2.TF>  0141*0 

CALL  F’PLOT  {  8 .  ZU  r  Z  .  1  »NSYP. TE  >  014170 

00  TO  100  01 4 1  BO 

C  WRITE  TEMPERATURE  OF  MEASURED  SPECTRA  0141*0 

1071  CONTINUE  014200 

CALL  SYMBOL  <  0  .  .8.3.  .  2.31  H  TEMPER  A  TUBE  -  LOR-  HEG  K  ,  0 .  .  .11  1  014210 

CALL  NUMBER <2. 6.8. 5*. 2. TF . 0 . . - 1 >  014220- 

CALL  NUMBER! 4.A»B»5».2.TE»0.»-t>  014230 

C  WRITE  SCAN  NUMBERS  014240 

YT=9,?9  014250 

XT*2 ,  014240 

CALL  SYMBOL <0.»9.29f.l9» 9HSCAN  NOS . . 0 . . 9 )  014270 

PO  1080  R-M.K2  014280 

FPN.IPLIST(K)  014290 

CALI  NUMBER  1XT.YT..19. FPN  *  0 . . - 1 )  0)4300 

XT-XTL1.  014310 

IF<XT.LT.12. ICO  TO  1080  014320 

XT«2.  014330 

YT-YT-.3S  014340 

1080  CONTINUE  014350 

C  LABLE  SLIT  SI7E  0143*0 

IF  < M ICRON. EO . 250 1C ALL  SYMBOL <8 . . 8.5. . 2 . 20MSL IT  •  2.5  ANGSTROMS. 0. .014370 

120)  014380 

IF(MICR0N.E0.251 1CALL  SYMBOL <8 . .8. 5. . 2.20HSL IT  *  3.S  ANGSTROMS .0. .0 1 4390 

120)  014400 

IF ( M ICRON. EO. 450) CALL  SYMBOL < 8.  .8 .5.  . 2 . 20HSLIT  -  4.0  ANGKTRI1MS . 0 .  . 01  44 10 

120)  014420 

100  RETURN  014430 

ENO  014440 

C  014460 

C  014470 

C  19  PRINT  PLOT  0I44P0 

C  014490 

C  014500 

SUBROUTINE  PPLOT (KOP.X. Y .N1.N2.TITL)  014510 

DIMENSION  TI TLE< 70) »X<!).Y(1)  014320 

DATA  HCTR/4840 . /. MGRD/12./. VCTR/ 4 ./ . VGRD/1 .71437 .FST/0.7  014530 

DATA  TITLE/6HWAUELN.AH  < ANG) . 6HAMPL I T . 3HUDE . - 1 . >0. .BM  FILE.  014340 

♦8M  SCAN. 8)<  SUMMED. 0.  .0.  .0.  .8HTEMTTRAT. 3HURE .0.  .0.  .  014350 

+8H!K>  EOR-.O. ,0. .O.ilHE.O. .0. .1X0.0, ,0 . . 1HC.0 . »0. , 1HP.  0145*0 

+0..0..1HA/  014570 

IE1FST.NE.0.  )G0  TO  11  0145110 

DO  10  J=35.70  014590 

10  TITL£!J>*0.  014*00 

DO  101  JbSO.62.3  014*10 

101  TITLE<J)--1.  014*20 

PO  102  J-40.42  014*30 

102  TTTLE<J>»-1.  014*40 

TITI  F148MTITI.E !  4*>»1  •  014*50 

1 1 n r  <64>-999.  014**0 

Tin  r<34)-3IIXXX  014*70 

11  FST=1,  014680 

K0M=3  014*90 

IF<Kt)D.FQ.l>X0N-0  014700 

IF!K0D.E0.8)K0M»4  014710 

TITLE  <  63 ) *K0B-1  014720 

IF<K0D.E0.7)T.TLE<*5)-29.  014730 

IFEKOD.NE , 1 >G0  TO  12  014740 

T I TLE !10) *N1  »  TITLE! 1 1 >*N2  *  TITLE! 121-TITL  01471.0 

12  CONTINUE  014760 

IFIR0D.E0.7)TITLE< 1*)«TITL  014770 

IF!K0D.E0.8)TITLE!18)-TITL  014  7110 

IF!K0D.GE.2.AND.K0D.LE.*) TITLE! -3*K0P43B)9TITL  0147*0 

CAU  MATRIX. Y.N1.N2,  ROM,  MCTR.HGRD.WCTR.  VGRD.  TITLE)  0141100 

RETURN  014010 

F NP  014020 
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SUBROUTINE  SUPAPL ( IV* TYPE . IAXIS* HUH* A, 10IM. ID) 

MAKES  VARIOUS  FORMS  OF  PLOTS  FOR  ALO  TEMPERATURE  PROBLEMS 
IV-O.l  FOR  DELTA  V  •  0  OR  DELTA  V  -  1 
T YPF  *0*2* 5.F0R  A  I  INEAR. 2  CYCLE  LOG.OR  SQUARE  ROOT  PLOTSiRESPEC TIVELY 
I AXIS-O* 1  FOR  NO  NEW  AXIS  OR  NEW  AXIS 
NUMPA=NO.  OF  PLOTS  TO  BE  MADE < NOT  NO  OF  AXIS) 

NUMBA  5  OR  LESS  IS  PEN  AND  INK  PLOT  .  6  OR  GREATER  IS  CRT  PL 

I  DIM-' DIMENSION  OF  INPUT  ARRAY ,  YARRAY .  FROM  LEFT  TO  RIGHT 
rr  IDIM-0  NO  CURVE  IB  PLOTTED 
YARRAY-ARRAY  TO  BE  PLOTED 

ID-3  IDENTIFICATION  NUMBER  TO  BE  PRZMTED  WITH  AXIS<4  DIGITS) 

IT  1 D=0  NO  NUMBER  IS  PRINTED  AT  ALL 

INTEGER  TYPE 
DIMENSION  Y ARRAY  <  t  > 

DIMENSION  PR  IDO) 

PR  ID ( 1 )  10HR 1 1 ROSSER 
PR I D ( ? ) «DA  TE  <  DDUM ) 

PRID<J)«T IMC(TDUM) 

IF < MSI AR I. NF. 362432 1041. AND. NUMBA. LE. 5) CALL  PLTID3(PRID. 340. * 1 1 . , 

!  1 .  ) 

IF  (MS  I  ART  .NE  .  3424321061  •  AND.  NUMBA  .LE  . 5)CALL  NEWPENO) 

IF (MSTART  «NE . 3624321061 . AND. NUMB A. GE .6)CALL  CRTPLT CPRID. 1 . . 1 7 . > 

IFCMSTART.NE .3624321061 >HEND*0 

IF (MSTART. NE. 362432 10A1 )CALL  PLOT ( 1 . * i . . -3) 

MSTART* 362432 l 061 
MEND-MENPU 

IF<IAXIS.EQ.0)GQ  TO  400 
C  AXIS 

IF(MEND.NE.1.AND.NUMBA.GE.6>CALL  FRAME (I. 5*1 .5) 

IF (MEND.NE. 1 .AND .NUMBA. LE »5>CALL  PL0T<17. *0. *-3> 

C  THE  FOUR  BORDER  AXIS 

CALL  PI 0T(0. *0. *3) 

CALL  ( IOT<12.»0.*2) 

CALL  PI  0TU2.  *6.  *2) 

CALL  r*L0T(0.  *8.  *2) 

CALL  Pl.OT (0.  *0.  *2) 

L  X  AXIS  FEDUCIARIES 

DO  120  K-l*2 
DO  120  1-1*149 
A-  *  05 

IF< I .E0.34< 1/5) )A- . 1 
IF ( I . EQ. 104(1/10) >A-. 13 
P-I 

P-.O04B 

C-0. 

IF<K . E0.2>A-9. -A 
IF(K. tO. 2)C-0. 

CALL  PLOT (B* A*3) 

CALL  riOT(B*C.2) 

120  CONTI NUf 
(  LADLE  X  AXIS 

IA- 4HOO-2004IV 
I B  •  I A I  ISO 
DO  303  M* IA* IB* 10 
X  M-  I A 

IF (M. LO. 504 (M/50) )G0  TO  301 
FPN*M  I004(M/100) 

X*. 0R4X-.il 

CALL  NUMBER (X.-. 35* . 13.FPN.0. • -1 ) 

GO  TO  30? 

.101  FTN  M 

X~.0H4X-.35 

CALI  NllMBFR <  X  •  -  •  4 »  ,2*FPN»0.  *  1 ) 

302  CONTINUL 

303  CONTINUE 

CALL  SYMBOL <3.8* ’ .8. ,2*2?HUAVELENGrH  (ANGSTROMS) *0. *22 > 

CAI  L  GYMDOI  (6.2*-  .6* .2* 1H. *0. *  1 > 

CALL  SYMBOL  <7.4*-.0».2»lH/*O.*i> 

CAP  I  P‘1.0 1  <  3 .  H/f  1 . 25*. 1) 

CAU  1101(4.07*  1.25*2) 

CALL  II  (It  <3.  V/*  1.07*2) 

CALL  PL0M3. 07.-1.25*2) 

CALI  SYMBOL <4. 12* -1.25*. 25. I4HV-  SEQUENCE  A*0.*14> 

CALL  SYMBOL ( 7 . 62 *  1 . 25 * . 15. 1ML . 0 . . 1 ) 

CALI.  SYMBOL  (  7.77.  -1.25*  .25*  lWO.O.  .1 ) 

FF'N-*  I V 

CAU  NUMBER! 4.62**1.25*. 25*FPN.0« .-!> 

C  l  API  I  Y  AXIS 

Tf  < IP)  1004.  1004# 1005 
1004  CONUNI/r 

CALL  SYMBCH  <  (.1*2.3* .2. 26HRELATIVE  INTENSITY  *90.*26> 


0JS3-C 
0153**' 
01  r-Z  * 

{  { •’  <  i> 

o  i  r  ■  4  i 
•>1V4  !  ' 
'U.K  o 
ol  •  ■>  • 
o  r .  4 
'»r-s 

*  i :  ' »  * 
DID 4  *'• 
0|'.-|:0 
OK- 1>*0 
M 

oii.r.p* 
<»ir-L.  . 

Aj;,* 

or-  •  • 
or.5» 

or/ 
or.:  • » 

■  i  s:  •  • ' 

or.*, 
or. 
OK-  U 
0  I  Sr- 20 
•MS.-30 
015640 

0156*0 
01  Sr  6  ■> 

ft  I  J,  N  "l  > 

or  v. 
015/  • 
i  ' 
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00  TO  1006 

1005  CONTINUE 
FIDNO-ID 

CALL  SYMBOL  l-laltl.li « 2» 26HRELAT IVE  INTENSITY  NUMBER » VO . r 76 ) 
CALL  NOMBERt-1 .1 ,6.3» .2»FlDN0r90. r-J ) 

1006  CONTINUE 

C  LINEAR  Y  AXIS  FEDUCIARIES 
IF (T*PE *2 >1020 *200. 1021 
1020  DO  180  K- 1  » 2 
DO  lflO  1-1,99 
A-.OS 

IF<1  .Ef>.5*<  1/5) >  A«  < 1 
IF<I .Ed. 10* < 1/10) >A-, 15 
B«I 

B-.O04B 

c*o. 

IF<K.E0.2>A*12.  A 
IF (K . EQ. 2)C'l2. 

CALL  PLOT l  A ,  P  r  3 ) 

CALL  PL0T<CrB»2) 
l 80  CONTINUE 

C  LABLE  Y  AXIS  linearly 

CALL  NUMBER <-*?5*-,l»,2*0# »0. ,  -1  > 

DO  340  M*10»90» 10 
Y*N 

YY-.08*Y-.l 

CALL  NUMBERS -,45,YY» . 2»  Y  »  0 . r-1 ) 

340  CONTINUE 

CALL  NUMBER;- .65» 7. 9f ,2» 100, ,0, ) 

GO  TO  270 

C  LOGARITHMIC  Y  AXIS  FEDUCIARIES 

200  DO  260  .1*1*2 
DO  260  K  * l , 2 
DO  210  MA-1,4 

DO  210  N«l»4 
ANA=MA- 1 
A*M 

A*4. *ALOGlO( . l*A4l .4.54AHA) 

FK--K 

A=-AE4  *  4FK-4 . 

FJ*J 

B~ 1 1 . 954FJ-I 1 ,95 
C=*P»  ,05 

CALL  FL0T<B.A,3> 

CALL  FL0T<C»A»2) 

210  CONTINUE 

DO  220  MA>3,6 
DO  270  M- 2,8» 2 
A*M 

AHA- HA 

A=4 » 4AL00101  •  1 4A4-AMA) 

FK*K 

A-AM.4FK-4. 

F  J~  J 

B*1 1 . 904F J-l 1 , 90 
C^B  f ♦ ! 

CALL  PLOT (B» Ar 3) 

CALL  PLOT (Ct Ar2) 

220  CONTINUE 

DO  230  M= 1 » 8 
A-M 

A*4  #  4A1.0G1 0  <  A  + 1  .  > 

FK-K 

A‘A44.*FN-4, 

FJ--J 

B*1 1 .0O4FJ-1 1 .80 
op*  .2 

CAll  PLOT  <  Pr  Ar.l  > 

CALL  PLOT <C  » A,2  > 

730  CONTINUE 

DO  740  M* 1 » 7 
A-M 

A* 4#  *AL0G10( A+ . 5) 

FK*K 

A*A  M . BFK-4 « 

FJ»J 

B*1 1 ,05*FJ-1 1 ,P5 
C*Bf . 15 

.CALL  PLOT < B* A« 3 ) 

CALL  FL0T(C*A,2) 

240  CONTINUE 

DO  250  M-M»3 

Ar  M 

A3  4  •  4  AL  OH  1 0  T  A  f  6  •  5 ) 


015/10 

015?.' 

otr.,-  io 

015  40 

015.'7  > 
•MS''. 

oi r-  - 

0  l  5 
•*>15  ■ 
i  f  ■ 

»'  I  '  ! 

<■  \  '  ■  '•« 
<M',h  ; 
01*.  t«. 

■  ML  *  ■ 
01  l-l  • 
<)  f  r.j  •  • 

cj  V  «.* 
o  J  *  •  •  VO 

or.  .* » 
'T'-  \> 
or. 

0 1 '»  ■“  '  ■ 
O I  V  l«‘ 
ojr.^  .o 
or.i'»-  » 

or.  -o 
015  •  • 

or  • 

01<V 
01  <0 !  •' 
0 !  f'1,  ■ 

O  J  A  ■  ■  * 
0  l  6'-4  ' 
0 !  t  «'• 
010  s  * 
0)6.  /•> 
016000 
Or.'  v. 

HM'Ti 
“1M  t- 
C  1  .M  . 
o  I  .*  ' 

0  16 

'  I  *•  I*  .* 
lis:  • 

<  1  M  *  ’ 
OlMl 
.1.1  ' 
OJ6. 

•M  ••  I 
t  *» 

•  I 

.*16  ’ 

)  s 


■  '  l  6 

or 
01/ 
o  j  /. 


o  i  •• » . 


0  !  •  •  '•  •  * 
O  |  ..  <> 

•  1  /  t 

•'  1  6  >  .0 
0  16  >  v*> 
| >  *r. 
01  VV.o 

o  r.  .  ’ 

o  i  /•  i 

o  1 •  '  ♦ 

O  1  6  •  ■ 

•  I  ' 

•M  '  1  • 

■  »| 

'I'  f 
01  4  ■ 


61 


Fk-k 

01*550 

A-A44.4FK-4. 

0165*0 

FJ»J 

016570 

b-ii  .estrj-u.BS 

0165BO 

C  =  R4 . 15 

016590 

CAt L  r LOTI*, A, 3) 

016600 

LAM  PLOT!!',  A,  2* 

0  1  66  1  «* 

r*»o 

CON  1  INUE 

016c . O 

A  4. 

016< 

F  »-J 

Ole'  4* 

B-  1 1  •/■J4FJ-1  i  » 75 

01 6  ‘  • 

C‘P4.?5 

Ot  6  -  -  ■ 

CALI.  PI  Of  !B*A.J> 

Of  6<*»  •> 

CALL  PI  OT  <  C  •  A » 2  > 

01  <.  .:M 

:»4o 

CONTINUE 

Ol  /./ 

l API C  V  AXIS  LOOARYTHMICAM  Y 

0  1  ^ 

PO  310  H - t , 9 

OlA^I  • 

f  I  N*  M 

0 1  6  7  ?  v 

Y  -  4 . 4AL0G1 0! FPN  > - , 08 

Ot*7x  ' 

CALL  NUMBER! -.2, Y, . 15»FPW,0. .-1 > 

01/ 

MO 

CONTINUE 

0167. 0 

f‘0  320  h*l 0,90, 10 

016 

FFN •  M 

01c  '  ) 

Y=4.*AL0G10<rPN,-.08 

**16  PO 

CAl  L  NUMPF. R <  -  » 35*  Y »  .  llWTI  N,0.  »-l) 

016' 

i:»o 

(UNI INUE 

oiau-v 

CAl  1.  NUMBl  Rf  -.5,7.93,  .  IS#  100.  *0.  »  -  l  > 

0 1  .‘>0 1  0 

no  TO  270 

01  AL  .'0 

10J1 

CONI INUE 

016??0 

1  API  F  Y  AXIS  FOR  SOUARE  ROOT 

0  1  6  6  J  0 

CAl  1  NUMBER! -•  2, -.08# .15,0.  ,0.  ,  1 ) 

016  '  >0 

PH  1022  M- 1 0 »  90  ,10 

01  6".  60 

f  r  N 

0160/0 

Y  . fl*SQRT  <  FFN ) - . 08 

OlAfif o 

CALL  NUMBER ( - • 33 » Y , . 15, FFN, 0 . » - 1 ) 

016G?  • 

1022 

CONTINUE 

01690-' 

CAM  NUMBER  < -,5,7.97, . 15* 100. ,0. ,-l> 

016910 

SOUARE  ROOT  FEPUCI ARIES 

01692C 

DO  1028  K»1 ,2 

016930 

BO  1028  1*1,99 

016940 

A  - .  05 

016^50 

IF  (MOD!  1,5) .EO.OlA-.l 

01696' 

I F <  MOD< 1 ,10) ,£Q,0>A~ . 15 

0t69-0 

hi 

0169-.  ' 

B-.fl*SORT!B> 

016'*  • 

r  -  o. 

01 7001 

If <  K , t 0.2)A-12. -A 

017010 

IF  (K.EO.  2)012. 

0170?  ' 

TALI  PLOT ( A , B , 3 ) 

01  7.  V) 

CAM  PLOT  ( C ,  B,  2 ) 

01 70 4 v 

io?a 

CONI INUE 

0170*  • 

1*70 

CONTINUE 

01  •  6U 

400 

CONT INUE 

IF  (  I fiFM.E 0.0)00  TO  407 

('1  V!  .1 

PI  151 

oi  ' 

17-3 

’I  ' 

Y MAX  -  0. 

oi  -»M  • 

no  401  !*1»  HUM 

■  1  -  ’ . 

FF! Y ARRAY! I ) .OF . YMAX 1 YMAX :  Y ARRAY! I) 

0  I  /  1  •  • 

401 

CONTINUE 

017  14 

DO  405  J  -  1  ■  I P I M 

V  1  t 

"  1  7  1  ’  ' 

..  J  ' 

*  1  1 

DIM  1  DIM-1 

1.5  • 

X  ;  »  4  <  17. /DIM) 

0|"l  ‘ 

If  (  1  FI  F  2  >  407  403,40/, 

Cl  - 1  •  • 

40? 

Y* YARRAY  (  I  )*<(). /YMAX  > 

01  ’  •  O 

GO  IIJ  404 

01  '  1  0 

4  03 

JF!  r ARRAY! I ).l  1  .YHAX/100. ) YARRnY! I )  -  YMAX/ 100. 

o  j  7:\  • 

>  4  ,  tALOG  1 0!  100.4!  YARRAY !  1  >/YM*,X>) 

■  i  . ' 

no  rn  404 

M  .  >- 

406 

IT (YARRAY! I ) .{  T.O. >  YARRAY (  I )  *0 . 

.‘t  • 

Y  *6 . *SOR I (YARRAY ( I > / YMAX > 

01 

404 

IF  <  Y  .  If  .  0 .  >Y=0, 

Ot  »' 

CAl.  1.  I  I  01  (X,  Y,  IT> 

01  7.T  * 

IT? 

01  - '  * 

4or. 

CONI INUE 

•  ‘5  ■ 

407 

RF I URN 

0  |  » 

END 

f»l/ 
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SUBROUT INF  APLQT <NEiND ) 

DIMENSION  IF0<3).B<12I> 

COMMON  Y f 300 >  *X( 300 ) *  YW(300) *  Z<  001 >»ZU(801>#T<19)» DUH<  801 ) 
DATA  1T0/10H( IX* 13* IX* » 10HF3.2* 1X*12* 10M1A1  )/ 

NE 'NUMBER  Of  ELEMENTS  IN  ARRAY  TO  Pfc  PLOTTED. 

Y “ARRAY  TO  BF  PLOTTED. 

ND-NllHRffi  OF  Ilf  CARES  FOR  LOO  RIOT. 

ND^O  HAMS  I  INFAR  PLOT. 

ND-Nf.'G  A HUE  GIVES  THAT  FRACTION  OF  A  Hf  C AW  .  . 

CALL  Ml  OT<-l,ltUHHY»NSIZF)  TO  SET  WIDTH  Of  PLOT  • 

DEFAULT  VALUE  IS  120  SPACFS. 

MAXIMUM  urcru  IS  120  SPACES. 

LATA  A 0/1H  /.Al/lH./ 

NSI7E  20 
SIZE*NS17f 
YM  a0 . 

PO  10  N= l »NE 
YN  -ANAX 1(YH*Y<N)) 

IF<NP.NE.0)L*N5IZE/ND 
lFLNP.l  T.O)F »-NSI2E*ND 
DO  40  N»l*N£ 

YN=Y(N)/YH 
IF (NP.EQ. 0)G0  TO  21 
YN-A MAXI <  YN* . 000000001 > 

NY  - AL0G10T  YN>*F4SIZE41 .000001 
GO  TO  22 

NY  SIZE4YNL1 .000001 
NY=  NAXO< 1 *NY ) 

NY*MIN0<NSIZE<1 *NY) 

NSIZEI-NSIZEU 
PO  30  J«1*NSIZE1 
D<J>*AO 
B<NY )-Al 
YN=ABS<  YIN ) ) 

YN--AMAX1  (  YN*  .000000000001 ) 

IAL*ALOG10<YN)F1.00000000ai 
IALr  31 ~IAL 
I  At. -HI  NO  <  31 » IAL  ) 

1 F  < I AL  *  LE • 26  >  GO  TO  37 
CALL  HXPUTX< IAL* IFG< 2 )  *  4 ) 

WRI  TF.  <  A*  IFO)N*Y(N)#(B(J>*J«t,NY> 

GO  TO  40 
NYN  Y(N> 

UR  1  IE ( 6*  292)N*NYN* J*J » NY ) 

’  FORMAT  <1X»I3.1X*I5*1X*121A1> 

CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONT INUE 
CONTINUE 
RETURN 
END 


SUBROUTINE  AL2<N1 *N2»NSTR) 

DIMENSION  NSTR(B) 

COMMON  Y(3OO)*XT3OO)*rU13OO)*Z(0Ol)»?U(0Ol  >*T<19)  *DUM<0O1 ) 
URI TE  <6*  201 ) 

5 1  rORMA  T ( 1 X • 1 3HENTER  PROGRAM  > 

NUN-  1 

DO  110  HN^NlwN 2 

CALI  ISNN(NN*  41  *  I DUMMY) 

DO  100  N^l*  149 

I«300*N4N5TR<NN>-301 

JSNNAX~0 

DO  20  K*1 • 300 

II=I4K 

CALL  ISNN<ZI<<1  TSN) 

IF  < I SNMAX *  L  T , 1 SN  > 1 SNMAX  ^ I SN 
20  Y(K  >* ISN 

IF (ISNMAX.LF.O)OO  TO  105 
CALL  SMOOTH <  300*11) 

NIP* 10006NN4N 
CAtLOUrU  NljM  »NID*  JOO ) 

ioo  continue 

105  CONTINUl 

PRINT  3<>I*NN*N*NHM 


017360 

017J70 

0173PO 

017390 

017400 

017410 

01  /A?0 

Ol  ?4«* 

01 /44A 

01  74?. 

01 7400 

0  l  74 70 

0 1  741*0 

0134V" 

017500 

01  7510 

017520 

017530 

01 7540 

0  I  7550 

017560 

017570 

0 1  7580 

01 7590 

01 7600 

01 761 0 

01 7620 

017630 

017640 

017650 

017660 

017670 

0 J  7600 

017690 

017700 

017710 

017720 

017730 

017740 

017750 

017760 

01 7770 

017700 

01 7790 

017000 

Ol 7010 

o l 7030 

017030 

017040 

017850 

017860 

01 7P70 

01 7000 

01 7090 

01 7900 

017910 

017930 

017940 

017950 

01 7960 

017970 

017900 

Ol  /v»9f> 

oinooo 
0  1 110 1  0 
0111020 
0100.10 
010040 

OUPV.0 
0 l HO 60 
0100/0 
01 0000 
0 l 0090 
010100 
018110 
010120 
018130 

0101 40 

0111150 
018160 
0181 70 
0181 0A 
0101  •>') 
01 11300 


! 


•4*  ft 


201  TOKHATl 1X.3IA) 

1 10  CONTINUE 

CAiL  ENDPLT 
RETURN 
ENP 
C 
C 

C  ?3  SMOOTH 

i: 

t: 

SURROUUNE  SMOOTH!  NDIM»NSN> 

COMMON  Y< 300>*X!300>* YU<J00VZ<801M7W<80I )»T(l9)rS<e01 > 

IF<N«H.lf .!>  00  TO  100 

N^NSM/2 

FNSM  *N«2+1 

DO  10  K  1 »N 

K2-N1R4 I 

6<M1  >«Y!KJ) 
to  S!K?>3Y<N> 

DO  AO  !=■  1  *NP  IM 
DO  20  K-1*N 
20  S<K>“S!K  H) 

s<nm  >«y< i> 

Y <  I  > -0 . 

DO  30  K~1  *N 

IKIfK.GE.NDrMM )  00  TO  40 
30  Y( I )*V< I >4Y< I4K) 

40  N1 =Nf  I  NDIMM 
N I -  MAXO <  N 1  *  I ) 

N1  ~NfNt 
DO  50  K3 1  »  N 1 
50  Y<  I  >  =  Y<  I )  ISO*  > 

AO  Y  < I ) =Y<  I ) /FNSM 
100  RETURN 
END 
C 
C 

C  74  OUf-L 

i: 

i; 

SIIPRfKJI  INC  QUPL  <  NO* ID*  ZE  ) 

I!  Y  ARRAY  IN 

I  NO  OF  PLOTS  OUTrUTTED  EXCEPT  MAKi  NO-1  AT  FIRST  ENTRY 
(  ID  OF  nor  IN 
C  r  7E  OF  Y  ARRAY 

DIMENSION  PROGXP! 3> 

COMMON  Y ( 300) *  X  < 300 )  *  YU! 300 ) *  Z ( 001 ) *  ZU( 80 1 )  *  T  <  1 9  )  * DOM ( 801 ) 
F  ROG 1 D ( 1 >  - 1 OHK I TROSSER 
PROG  T  D  <  2 ) *  DA  TE ( DDUM ) 

PROG  I D  <  3  >  »  T I  ME  <  TDUM ) 

iriNo.ro. i >call  crtplt <prooip*i . *  t7. > 

irtNO.FO.l >CAtL  PLOT! 3 . 45  *  0 . * -3 ) 

I F  <  NO •  F  0  •  1  ) GO  TO  At 

I F ( NO  !NO/4>*4.E(1. 1  >CALl  FRAME! 4. 15* .5) 

IF ( NO- < NO/4 ) 44 . EQ« 1 ) GO  TO  A1 
CALL  Pl.OT  (4*25*0.  *  -3) 

A 1  FID::  IP 

CALL  PLOT <  -4. , .4*3) 

CAI.I  PLOT (0. * .4,2) 

CALI  PI  OT  ( 0 ,  *  1 0. 4  •  2  ) 

YMAX-0 • 

no  j  i »  t^e 

IF(Y  <I).QT.YMAX ) YMAX’Y  (J) 

l  CONTINUE 

rr  <  YHAX.I  f  .0.  )G0  TO  3 

CAI  »  NUMPER!  2. *-.4, .2, FID. 0. *1 > 

CALL  SYMDOl (-3. *0 . * . 2* 3MYMAX- *0 . *5) 

CAM.  NIIHDER!-?  ,  *0.  *  «  2*  YHAX  *0.  *  I> 

113 

DO  7  I *! *  ?7E 
V Y ”  4  »  4 Y  <I)/YMAX 

ir< yy.gt .o. >  yy • o. 

FI  -*1  —  1 
FI7  IZF  I 
<x=.44io.*r i/r i z 
(  ALL  PI  01  <  YYrXX,  IT  > 

II  *7 

CUNT  TNMI 
NO  NO  I  I 

3  LON  I  I  NUT 

Rf 1  URN 
(  NIi 


018210 
018270 
018230 
018240 
018250 
018270 
010200 
01 82 VO 
010300 
010310 
010370 
018330 
016340 
010350 
010340 
0103/0 
010300 
018390 
010400 
010410 
010470 
018430 
018440 
018450 
016440 
010470 
018480 
018490 
01G500 
Cl  O'.  10 
010520 
018530 
01 6540 
010557 
0181  AO 
010570 
010590 
018400 
01 84  JO 
018470 
010430 
018440 
010450 
018440 
01 04 70 
018480 
01  PA’.O 
0J8700 
010-10 
018720 
018730 
010740 
010750 
016760 
01 8770 
019780 
018*90 
018900 
01 00 10 
018070 
0100*0 
010840 
016050 
010860 
010070 
010080 
01009' 
018900 
018910 
010970 
018930 
Ot  0940 
010950 
018960 
010970 
018900 
018990 
019000 
01 9010 
019070 
0 1  90  V) 
0190*10 
0190V' 
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L  0193*'’ 

r  ?A  SI  IT  PfcllfSKAtt  OJV.r.’.' 

r  O I ?.W 

f;  oj  vt vo 

PROGRAM  SI  I f< INPUT, OUTPUT .TArt/e  I  Af-*f  A  » Oil  I  Pll  T  »  T  APF  7  ) 

IUMFNS10N  X(l00>»2<100>.V<l0>'|iM0>’ii<  JO)'PA(.l>  oi  ^io 

PRINT  909  0|V4;o 

VO?  FORMAff/*  •••••RPHTHIH  k  TO  rOf’V  tAPf  7  AFTER  EXECOMt  'MvUc 

*  *0N.<«*§*  */>  Of,  J.U 

RE u ! Ni»  2  o i v 

RF Alt T  L*  •  V02  I  AAN  0 1  v .!/■'> 

NAANf  .000 >01 

90."'  F0RMATTF9.3 >  0 1  ■ 

Ml  K;  l  at  ,N  0  1  v*90 


<  1 1 1  ^  I 

•/IV.'"' 

KFAn<2«902MZ<n»l'»t»N> 

0  1  VS |0 

URI  IfTAf 904) 

0|7V« 

V  *4 

F  ORHAT  <  // JX*X* 3X*Z*/  > 

01  v.  O' 

UKl  re.<A#  V03)  <X<  l  >  »/<  T  >f  1  l  *  N ) 

o  |  • 

VO  1 

1  l)R MAT  <  IX»*  J.Op  IXfF'5.0) 

'•1  ,  •'» 

r 

UT  1  ►  I  I'ri’F-C  0(3)  U  U<4>A 

0  1 

V'.»>*  .54T/T ! ) f?(N> > 

0* 

V<4 1*0. 

oi  v:,». 

Ml  10  1-1 ,N 

oj  vr.vu 

It  <Z(J)  .Gf.V(4>  FJXaJ 

01 «AOO 

IF  <  l <  T  > .  fi  r . 0 <  * )  > V<  *  >  7 T t ) 

0 1  9  a  1  •) 

to 

rUNTTNt  IE 

oi9a:*o 

UT  4)--V<4)  -v<n> 

0  |  9A  ?0 

VC 1 >*X< JX* 

A19A4.'- 

7U  VC)H  .  :,*0<4) 

0  1 9  .•  *•••  * 

Ml  IS  1  l-N 

OJ  9A-0 

JX»  I 

0 l 9A  70 

II  (  /W.l  T . /(  I  )  ) GO  ro  1A 

0194P  1 

r< 

(UN  f  INI  IE 

019*90 

1  A 

mi  I?  t».m,N 

0  |  9  -On 

jx  ?-  r 

0 1  c  >  |  A 
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